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Abstract 

Linear space-time block codes (STBCs) of unitary rate and full diversity, systematically constructed over arbitrary 
constellations for any number of transmit antennas are introduced. The codes are obtained by generalizing the existing 
ABBA STBCs, a.k.a quasi-orlhogonal STBCs (QO-STBCs). Furthermore, a fully orthogonal (symbol-by-symbol) 
decoder for the new generalized ABBA (GABBA) codes is provided. This remarkably low-complexity decoder relies 
on partition orthogonality properties of the code structure to decompose the received signal vector into lower-dimension 
tuples, each dependent only on certain subsets of the transmitted symbols. Orthogonal decodability results from the 
nested application of this technique, with no matrix inversion or iterative signal processing required. The exact bit- 
error-rate probability of GABBA codes over generalized fading channels with maximum likelihood (ML) decoding is 
evaluated analytically and compared against simulation results obtained with the proposed orthogonal decoder. The 
comparison reveals that the proposed GABBA solution, despite its very low complexity, achieves nearly the same 
performance of the bound corresponding to the ML-decoded system, especially in systems with large numbers of 
antennas. 

L Introduction 

The design of space-time codes (STCs) as building blocks of multiple-input-multiple-output (MIMO) systems 
[1] has been a popular research topic in wireless communications since the discovery of the Alamouti scheme [2] 
and the orthogonal space-time block-codes (OSTBC) [3], [4]. 

Although STC techniques were originally envisioned for rather modest spatial-temporal dimensions, the last few 
years has seen a growing interest on techniques suitable for large/scalable MIMO systems [5], [6]. Arguably, this 
recent owes mainly to two related developments. The first is the advance the antenna technology enabling the design 
of compact multi-mode antennas [7]-[9] which, though small in size, can provide a large number of independent 
channels by exploiting spatial, pattern and polarization diversities [10]. The second is the increasing popularity of 
distributed and cooperative systems [11]-[13]. 

These recent developments provide a rationale for the feasibility of MIMO systems of not only large but, in 
fact, variable dimensions. For example, in the case of compact multi-mode antennas, MIMO setups of different 
dimensions may occur because the number of channels effectively made available by the spatial, pattern and 
polarization diversities dependents on the physical scattering properties of the environment. Analogously, in the 
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case of distributed architectures, variable MIMO setups arise naturally since the cooperative transmission can only 
be employed after the group of cooperating devices is assembled [14]. 

MIMO techniques can significantly improve both user and rate capacities (efficiencies) of wireless communica- 
tions systems [15], [16]. Although space-time transmission diversity techniques in point-to-point links are known 
to be sub-optimum from a Shannon capacity viewpoint [17], these techniques have their extremely low-complexity 
to their advantage. This trade-off argument is especially compelling when considered that the capacity achieved 
by MIMO systems is in practice significantly limited by the use of digital modulation [18]. It is in distributed, 
cooperative diversity systems, however, that the flexibihty of scalable space-time coding techniques have the most 
of their potential [19], [20]. 

Given the above, it is fair to say that, when considering space-time block codes as building blocks for scalable 
MIMO systems, modulation-independence, generalized and systematic constructibility, orthogonal (symbol-by- 
symbol) decodability and rate-unitariness' are the most essential features to be sought. These requirements make 
a strong case for linear STBCs and since existing techniques such as the Alamouti scheme, the OSTBCs and the 
maximum-rate minimum-delay (MRMD) STBCs are either limited to a small number of antennas, or to particular 
constellations, or undermined by sub-unitary rateness [21], [22], many have taken upon the challenge of designing 
alternative solutions. To cite a few examples, the universal STCs by Gamal et al. [23], the full-rate full-diversity 
(FRFD) STCs by Xiaoli et al. [24], the full-rate fuU-diversity orthogonal STBC scheme by He et al. [25] and the 
variable-rate STBCs by Kim et al. [26]. 

Each of the above-mentioned methods are clever techniques representing a step forward in the quest for high- 
performing, low-complexity STBCs with no or little rate-loss. To the best of our knowledge, however, these and 
other similar contributions fail to simultaneously deliver a systematic design, scalable to any number of transmit 
antennas, applicable to arbitrary constellations and allowing a linear, symbol-by-symbol decoding method with 
performance comparable to that of the maximum likelihood (ML) decoder. 

In this paper, a family of space-time block codes that satisfies all these objectives is introduced. The codes are 
obtained by generalizing the ABBA (a.k.a ^Mosi-orthogonal) STBCs (QO-STBCs) presented in [27]-[29]. A linear, 
orthogonal decoder for the proposed generahzed ABBA (GABBA) STBCs is also provided. The performance of 
the new codes is also studied both analytically and through computer simulations. In particular, the exact bit error 
rate (BER) probability of GABBA codes over uncorrected memoryless fading channels with unequal powers and 
arbitrary statistics is derived under the assumption of perfect channel knowledge at the receiver and ML decoding. 
Comparing those results against computer simulations, it is shown that orthogonally decoded GABBA STBCs 
achieve nearly the same performance of achievable with the ML decoder (which is prohibitively complex for 
large codes and higher-order constellations), especially in systems with a relatively large number of antennas at the 
receiver. The transmission rates theoretically achieved by the GABBA codes are compared against the corresponding 

'Rate is here understood as the ratio between the number of transmissions epochs and the number of symbols transmitted per block. For a 
formal definition, please refer tolsl 
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Shannon SIMO bound [16], [30], showing that GABBA codes have the potential to reach the capacity of the SIMO 
channel, despite the remarkably low-complexity involved in the technique. 

II. Systematic Construction of GABBA-STBCs 

A. Preliminaries 

To the best of our knowledge, a simple and systematic method to construct linear, full(unitary)-rate', modulation- 
independent STBCs for any number of antennas is yet to be demonstrated. In [27], Tirkkonen et al. conjectured 
that larger ABBA STBCs codes could ''probably be constructed by repeatedly iterating the ABBA scheme.'" It will 
be shown here that this conjecture was almost correct, except that the generalized construction involves two distinct 
ABBA schemes. 

Consider a vector sk of symbols taken from a non-ambiguous arbitrary complex constellation S E C^^qj: 

SA'-[si,--- (1) 

where ^ denotes transpose. 

Definition 1 (Linear Space-time Block Encoding Matrix): 

The term linear space-time block encoding matrix will be used throughout this paper in reference to matrices of 
the form 

Cll • • ■ Ci,K 



eC^fo'jS (2) 



ck,i ■ ■ ■ Ck,K 

where each element relates to one and only one symbol sj. e s^, such that c^.j € {s^, — Sfc, s^., — s^.}. □ 
Notice that definition [Qrules out constructions based on linear combinations of entries of the symbol vector, such 

as the rate 3/4 code that appears in [3]. 

Definition 2 (Raw Symbol, Raw Encoding Matrix): 

The symbol Sk from which an entry c^j of a linear space-time block encoding matrix C{sk) derives, is referred 
to as the raw symbol of Cij, and is denoted by (ci.j). Likewise, the matrix {C{sk)) consisting of the raw symbols 
of the space-time block encoding matrix C{sk) is referred to as the raw encoding matrix of C(sa'). □ 

Definition 3 (Dense STBCs): 

A linear space-time block code matrix C(sa') is said to be dense if and only if (iff) (ci.j) 7^ V (i, j). □ 
Definition 4 (Complete STBCs): 

A linear space-time block code matrix C{sk) is said to be complete iff it satisfies both the following conditions: 

«) (c»,i(sA)) 7^ (Cj,i(sA)) V i^j, 

ii) (cx,,(sa)) ^ (ci j(s7<)) V i 7^ j . □ 
In plain words, if the rows and columns of an encoding channel matrix represent different transmit epochs 
and different transmit spatial channels, respectively, definition |3] implies that at every transmit epoch, all diversity 
branches are utilized. On the other hand, according to definition |3 the raw columns and raw rows of a complete 
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encoding matrix C{sk) are, simultaneously, mutually exclusive permutations of sk and s^, respectively. These 
two properties ensure that the code fully exploits the "raw diversity" of the channel, as defined in [27, section II, 
§1]. Notice that the K-hy-rit rectangular matrices resulting from column partitions of a dense and complete linear 
space-time block encoding matrix C{sk) preserve this property. For simplicity, such matrices will also be, therefore, 
referred to as dense and complete. 
Definition 5 (Rate): 

The rate of a linear space-time block code described by the matrix C{sk) is defined as the ratio between the 
number of independent raw symbols and the number of linearly independent rows in C(s;f). □ 

Notice that, according to definition [S] a linear space-time block code as defined here can, at most, achieve a 
unitary rate. Therefore, the terms rate-one and full-rate will be hereafter used interchangeably in reference to a 
Unear STBCs of rate 1. 

Definition 6 (Real Orthogonality): 

A pair of vectors a and b e C"^^ is said to be real orthogonal to one another iff 

n 

a'^-b = bT-a==^a,&, =0. (3) 

Definition|6lis introduced in order make clear the distinction between the orthogonality property given in equation 
Q and the notion of vector orthogonality which is usually defined, in complex Hilbert spaces, over the product of 
a vector against the transpose-conjugate of another [31]. 

Definition 7 (Block-orthogonal STBCs): 

A matrix C G C"^™ is said to be block-orthogonal iff there exists a permutation C of C such that 

Di 
D2 



(4) 



where ^ denotes Hermitian. 

Furthermore, C is said to be real block-orthogonal iff a permutation C exists such that 

Di 



D, 



(5) 

□ 



Definitions 6 and 7 establish a slight generalization of the concept of ^MOii-orthogonality, in the sense presented in 
[27]-[29], with q'Maii-orthogonality resulting in the special case where the dimensions of Di and D2 are identical. 
The notion of real orthogonality will be used later in the derivation of the GABBA orthogonal decoder, in section 

Definition 8 (Manifold Matrices): 

Consider the generating function defined below over the space of complex-valued square matrices A, B e C^^^, 



X(A,B) 



A B 
B A 



(6) 
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The 1'*' order ABBA manifold matrix of the vector S2 — [si, S2], denoted by X[gj g^]' defined as the outcome of 
S2), with Si e S. Accordingly, the 2"'^ order ABBA manifold matrix of the symbol vector 84 = [si, S2, S3, S4] 
with generating function X is defined as the outcome of ^(Xj^j S2]i^[s3.s4]) and denoted by X^^^ ^^j. 
In general, the n-th order ABBA manifold matrix of the vector s^^ with generating function X is defined as the 



outcome of X(X, 



(n-l) 



) and denoted by X^ , where Sx/2 is the complementary symbol given below 

T 



— [Sk/2+1-: ■ ■ 



(7) 

□ 



Lemma 1 (Commutativeness of n-th Order Manifold Matrices): 

The n-th order manifold matrices of independent, equal-sized vectors bk and b/f e C^^^ commute with respect 
to matrix product. 

Proof: Let us rewrite X^^ and XJ^^ in terms of their half block-partitions: 



-Xi"-'' X; 



1) 

bK/2 



DA72 

-x.^"-') x^ 



Ok/2. 





Ai 


Bi 




Bi 


Ai_ 


1 


A2 


B2 




-B2 


A2_ 



(8) 



(9) 



Writing the products X^^. 
obtain. 



XjJ^ and XJJ^ • X^^, explicitly and using standard block-partition properties we 



Ai A2 Bi B2 Ai B2 + B1 A2 
-(Ai-B2 + Bi-A2) Ai-A2-Bi-B2 



(10) 



X" -X" ^ 

Ok slk 



(11) 



A2-Ai-B2-Bi A2-Bi+B2-Ai 
-(A2-Bi + B2-Ai) A2-Ai-B2-Bi 
It is clear from equations dlOt and that the condition for the n-th order manifold matrices X!^^ and Xj^^^ to 
commute is that their minors Ai, Bi, A2 and B2 commute. These minors are, in turn, (n— l)-th order manifold 
matrices (see equations (|8ji and (|9}) and, therefore, commute under the same condition. Since the minors of first- 
order ABBA manifold matrices are complex scalars, the commutativeness of X^^ and Xj^^ with respect to matrix 
product follows immediately, by induction. ■ 



B. GABBA Encoding Matrices 

We can now introduce the major result of this section. 
Theorem 1 (Construction of GABBA Matrices): 

For any given number of transmit antennas nj e N+, there exists a dense, complete and block-orthogonal GABBA 
matrix for rate-one space-time block encoded transmission of symbols taken from any complex-valued arbitrary 
constellation S. 

Proof: The proof is constructive. For any given n^, let n — \\0g2 "^tl and K — 2", where \x\ denotes the 
smallest integer greater or equal to x. Next, construct the (n— l)-th order ABBA manifold matrix xi]J~^^ of the 
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symbol vector ™d the (n— l)-th ABBA order manifold matrix X^^^^^"* of the complimentary symbol vector 
both with the generating function X given in equation (|6j. 
Next, construct the dense and complete GABBA mother encoding matrix C{sk) by computing 



c(sK) = Yj;,=r(x(-i),xJ;f), 



(12) 



y(A,B)^ 



(13) 



with F(A, B) defined as 

A B 

-B" A" 

For obvious reasons, ^(A, B) is referred to as the wrapping function. 

Notice that the manifold matrices X*^"^^ are constructed over 2*^"^^ vectors each carrying a mutually exclusive 
partition of Sk tuples of symbols and, furthermore, that at each iteration of X, these arguments are in opposite 
quadrants of the resulting matrix. The same occurs with the final "wrapping" operation, i.e., the application of 
equation ( I13> . It is therefore clear that the GABBA mother encoding matrices are dense and complete. 

It is also evident that Yg"^ are square matrices of size K = 2". Therefore, the capital scalar K will, hereafter, 
exclusively denote integer powers of two. Rectangular (K-hy-nt) GABBA encoding matrix are obtained by simply 
selecting any nt columns of C{sk)- Since columnwise partition does not affect the density and completeness of 
C{sk), GABBA encoding matrices Q^^ are also dense and complete, which concludes the proof. ■ 

Hereafter we shall simplify the notation slightly and write Ca', instead of C{sk), for the K-hy-K GABBA 
mother encoding matrix. Accordingly, the matrix obtained by selecting n* columns of Ck will be denoted by C„t ■ 

The following property is an immediate consequence of Theorem ^ 

Corollary 1 (Block-orthogonality of GABBA Encoding Matrices): 

GABBA mother encoding matrices are quasi-orthogonal and GABBA encoding matrices are block-orthogonal. 
Proof: The K-hy-K GABBA mother encoding matrix is given by 



X 



X 



(n-1) 

SA72 





A 


B 




-B" 


A" 



(14) 



C(sa)-C(sk)" = 



(15) 



Writing the product C{sk)-C{sk)^ explicitly, and after trivial algebraic manipulations, we obtain 

AA" + BB" BA AB 
(B-A-A-Bf A-A^ + B-B" ' 
Equation il5\ indicates that the q'Mflii-orthogonality of GABBA mother matrices is conditioned on the commu- 
tativeness of its minors (A,B) and, therefore, follows directly from Lemma 1. 

Rectangular (K-hy-nt) GABBA matrices Qj^ are obtained by arbitrarily selecting rit columns of C(sa:), that is, 
by selection a total of nt column vectors out of the sets {c|_i(sa-), • • • , c^x/2{sk)} and {c^ i^xfJ,^K), ■ ■ ■ , ci^k{sk)}- 
Since this columnwise partition does not affect the orthogonality between these two sets, the corresponding mother 
encoding matrix, Qj^ are obviously fe/ocA:-orthogonal, which concludes the proof. ■ 
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An example of a GABBA mother encoding matrix with K = 32 and a Matlab implementation of the encoder 
described in this section are given in Appendix |l] The reader is encouraged to inspect the example provided and to 
utiUze the program in order to generate GABBA matrices of various sizes and verify their block-orthogonality^. 

III. Quasi-Orthogonal Decoding GABBA Codes 

A. Preliminaries 

Consider the problem of decoding the received signal vector corresponding to the transmission of a GABBA 
matrix C„( through a memoryless block-fading channel h,ij in the presence of an additive noise vector n^- At a 
given receive antenna, such a received signal vector can be modeled as 

r = ■Kt+T^K = Ha- • sk + n^, (16) 

where sk is the augmented symbol vector, defined as 

= \si,- ■ ■ ■.sk-.sI,- ■ ■ ,s\Y , (17) 

where * denotes complex conjugate. 

Examples of the use of alternative representations of the linear STBC receive vector, similar to that given in 
equation (I16> . can be found in the literature. For instance, in [32], the idea was employed in order derive a simpler 
decoder for OSTBCs, and in [29], the transformation was used to construct the ^MOiZ-orthogonal decoder for the 
4-by-4 code there proposed. 

We have shown, however, that this technique can be expanded well beyond the relatively straightforward applica- 
tions suggested in [32] and [29]. In particular, we have shown in [33], [34] that a more elaborate utilization of these 
transformation can be exploited to construct a novel decoder for OSTBCs which enables its use in time-selective 
fading channels. In this section, we build on the concepts laid in [33], [34] in order to show that GABBA codes 
are indeed orthogonally decodable. 

B. Structure of GABBA Encoded Channel Matrices 

Equation M6\ indicates that Hx are constructed by simply re-arranging the corresponding channel entries 
appropriately. These matrices are, therefore, hereafter referred to as the GABBA encoded channel matrices. 

Begin by noticing that when the transmitter has a number rit < K of antennas, an rif column-partition of the 
corresponding K-by-K GABBA mother encoding matrix is utilized. Mathematically, this is equivalent to the case 
when K — fit columns of Ck and rows of h (corresponding to the "unused" antennas) assume zero values-^, leading 

^Instructions on how to generate a GABBA encoding matrix of symbolic variables are given as comments in the Matlab code provided. 
Note that the Symbolics Tool Box is required. To encode a numerical symbol vector s, simply call the function as usual giving s as input. To 
inspect the block orthogonality of the GABBA encoding matrices, first generate a symbolic GABBA mother encoding inatrix <C> and then 
type <C' *C>. 

'This is the reason why we inaintain the notation in equation <16t 
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to a sparse Hk, but leaving unaffect its real ^wasj-orthogonality properties. It is, therefore, sufficient to consider 
mother matrices Hk when studying the properties of GABBA encoded channel matrices H„, . 

Before we proceed, let us introduce some additional terminology . 

Definition 9 (Notation for Matrix Partitions): 

The symbol [AJ^'p denotes the left (P=L), right (P=R), upper (P=U) or lower (P=D) partition of the matrix A, 
containing its p leftmost, rightmost,uppermost and lowermost columns/rows, respectively. 
Definition 10 (Notation for Matrix Reorientation): 

The symbols A'^ and A* denote, respectively, the bottom-up and the left-to-right reorientations of a matrix A 
(a.k.a. flip upside-down and flip left-right operations). Trivial properties of these operations are: 



where A e C"^™ and b e C™^^ 

Next, observe that due to the structures of Ck and s^, can be written as 

Hi^i 
1 

K — , 

Hk, 

where H/^j and H^a are both K/2-hy-K matrices. 
The following result can then be stated and proved. 
Lemma 2 (Construction of GABBA Encoded Channel Matrices): 

The non-zero minors and H-k^ of the GABBA encoded channel matrix Hk are given by 



(Aty = A, 

At-b = (A-b)^ = Att-bt 



(18b) 



(18a) 




vector h^, respectively. 




h^^ extended with K—nt zero-entries 



4 [[h+]D^^/2 [h+r^l^] = [hK/2+1, • • • , 0, • • • , • • • , , /II, • • • , /Ife, • • • , 0, • • • , hK,2]\ (23) 

^ ^ ' 



last K/2 entries of first K/2 entries of 



with the corresponding generating functions H and H given by 



A 



A B 



ir(A,B) 



(24) 



B -A 



iJ(A,B)4 



A -B 



(25) 



B A 
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Proof: The proof is constructive. The noise-independent terms of the first -ftr/2 entries of r can be written as 







-S2 Si 





lU:-^lD:2 



hi 

h2 



D:(K-2) 



hi h2 
/l2 -hi 



,U:2]R:(^f-2) 



[H^jD:(K/2-2) 



Si 
S2 



ID: (if -2) 



H 



-Ki ■ Sif . 



(26) 



Equation J26> derives directly from the structure of Ck (see section III-B> . and reveals that is the upper 
half-partition of the n-th order manifold matrix of h^, with the generating function given by equation (I24> . 
Analogously, the noise-independent terms of the last K/2 entries of r (flipped upside-down), can be written as 







'k_i 




2 


2 ^ 


L:f-2 




-i-'k 




2 


^ 2 



"-K2 



-hi -h2 
h2 -hi 



hx hx 

— + 1 — +9 

2^^ 2^- 

hn 
T+2 T+i 



-hK 



■K2iR:K2 



■K2-S 



A2 



(27) 



From equation i21i it is found that Hj^^ has the following structure: 



Hif 2 — 


hK/ 2+1 
hK/2+2 


-hK/2+2 
hK/ 2+1 


U: 2 

[HK2]R:^f-2 

L: if/ 2-2 


-/ll 

-/i2 


h2 
-hi 


Pk2]r|^/2-2 







(28) 



Equation (I28> finally reveals that Hj^^ is the upper half -partition of the n-th order manifold matrix of h^, with 
the generating function given in equation M5\ . which concludes the proof. ■ 

Lemma|2]gives an algorithm to construct GABBA encoded channel matrices of any size, and reveals the similarity 
in the construction of the GABBA mother encoding, and encoded channel matrices. As an example, the minors VLki 
and H/f2 of the encoded channel matrix corresponding to C32 are shown in Appendix HII A Matlab implementation 
of the algorithm is also provided. The reader is again encouraged to inspect the example, as well as to utilize the 
program to generate GABBA encoded channel matrices of various sizes and verify the veracity of equation M6\ . 



C. Properties of GABBA Encoded Channel Matrices 

Let us now introduce the following property of GABBA encoded channel matrices. 
Lemma 3 (Quasi-orthogonality of GABBA Encoded Channel Matrices): 

The following property holds and is referred to as the quasi-orthogonality of GABBA Encoded Channel Matrices"*: 



Hi 
H2 



(29) 



"^This represents a slight abuse of the terminology introduced earlier, given that the property is defined over a sum of matrix products. We 
believe, however, that this does not contribute to any confusion and simplifies that language somewhat. The reason why such a product sum of 
HT"^ and is considered is to became clear soon, when the procedure for orthogonal decoding is introduced. 
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where Hi and H2 are K/2-hy-K/2 matrices. 

Proof: Let Aj and be manifold matrices of the sub-vectors [hj 



+ lU:K/2 



and [hj 



OjU:K/2 



such that: 



a) Ai and Bi are structured as described by the generating function given in equation MA\ : 

b) A2 and B2 are structured as described by the generating function given in equation ( I25t : 

c) Ai depend on the iir/2-tuple ha = [/ii, • • • , hx/2], while B^ depend on the A72-tuple = [/1K/2+I7 ■ 
From equations ( I24t and ( I25t we then have 

^ Hi P 



where 



Hi = lAipH 



Ho 



P" H2 



Bi|2 + AT.A^+BT.B 
B^P + AVAt+BT-B^, 



2) 



"■l -^1 
P=Pl P«, 

Pi ^ A«.Bi+AT.B^, 



(30) 

(31) 
(32) 
(33) 
(34) 



Form equation ( I33> it is seen hat the ^Mai/-orthogonality of the GABBA encoded channel matrices, as defined 
in equation i29\ . can result from two independent, and both sufficient, conditions: 

i) Pi = P^ ("weak" condition) 
ii) Pi = ("strong" condition) 

The commutativeness of manifold matrices established by LemmaQJimplies^ that Pi = P|. It is therefore evident 
that the weak condition is only satisfied if either the matrices A^ and B,; or, equivalently, Ck are/is real-valued^. 

The strong condition ii, on the other hand, is generally satisfied, as shown below. 

Let Ai=i/(Ai,i,Bi,i),Bi=i?(Ai,2,Bi,2),A2 = #(A2,i,B2,i),B2 = if(A2^2,B2,2), where Ai,, and Bi,, 
are build over the the first and second half partitions of h^, respectively, while A2 j and B2.j are build over the 
first and second half partitions of h^, respectively. Then, Pi can be written as: 



dependent on 



[hi, 



dependent on 



Pi 



-^1,1 ^1.1 



BH 
1.1 



-A" 



Ai,2 

B1.2 



Bi,2 

-Ai,2 



■■^2,1 -^24 

A^ B^ 
-^2,1 ■'-'2,1 



B2,2 


A2,2 




Pl,l 


Pl.2 


A* 


^2,2_ 




P2,l 


P2,2 



(35) 



dependent on 

[''K/if I. -- •'•jf/2l 



dependent on 



^Although the generating function considered in LemmalTI differs from those given by equations I24i and 1251 . the extension of the proof to 
these generating functions is trivial. 

'in the later case because then the linear combination of the received signal vector can be performed independently over the in-phase and 
quadrature components of the channel vector hnt , splitting the corresponding encoded channel matrix Hx into real and imaginary parts. 
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From equations ( I29l l. i30\ . i34\ and ( I35t . and due to the properties a, b and c listed above, it is evident that the 
sum-product structure of P propagates to P, to Pi, to its blocks Pi,j, and so fourth. For example Pi i is given by 



Pl,l = (A'J i-Ai,2+B^ i-Bi,2 + B24-B2 2-A2 1-A2 2)-(A'J i-Ai,2+B5'_1-Bi,2+B2.i-B2 2-A2_i-A. 
same generating function 



* \H 
2.2) ' 



= Aq^i- Ai,2 - Ai,i ■ Aq_2 + ■ ■ ■ = A" ■ Bi - Bi ■ AH + ■ ■ ■ (36) 

dependent on dependent on same structure as Pi 

[hi ,-■ ,hK/4l ['»Jf/4f I--- .''Jf/Sl 

Ultimately, this structure propagates to each of the scalar elements p^^^ of Pi, yielding 

K/2 

i=i 

where the set of pairs {(/iai, hbi),' ■ ■ , {hiK/i^ ^^'^ different for each {n, m)-th element of the matrix. 

Equation ( I37t . together with the fact that Hi and H2 are matrices of the same size, implies the ^Maii-diagonality 
of P, which concludes the proof. ■ 

The reader is invited to utilize the Matlab code provided in Appendix III- Al to generate GABBA channel encoded 
matrices of various sizes and verify its quasi-orthogonality property algebraically^. 

D. Quasi-orthogonal Decodability 

The q'Mflii-orthogonality of GABBA encoded channel matrices has the following consequence. 
Theorem 2 (GABBA Q\i3&\- Orthogonal Decodability): 

The receive vector vk corresponding to a GABBA STBC K-by-nt transmitted through a memoryless block- 
fading channel with additive noise can be linearly combined so as to produce the following pair of vectors tk^ and 
fifj, respectively dependent on the mutually exclusive _ft72-tuple of symbols, as shown below 

fifi ^ilK/2 ■ SK/2 + ^K/2i , (38) 
^K2=H.K/2 ■ Sk/2 + nK/22 ■ (39) 

Proof: Consider the following linear combination: 

f^ = [H«f^.rK+rJ,.[H^f^. (40) 

^This is accomplished as follows. First, give the variable K an appropriate (power-of-two) value. For example, type <K=32;> in the 
command line (the brackets <> not included) to make K = 32. Next, generate the conesponding symbolic GABBA encoded channel 
matrix H/^ and its non-zero minors H^i and ^i'^ the command < [H, HI, H2 ] =GABBAEncodedChannelMatr ix ( [ ] , K) ; >. If 

desired, omit the semicolon or type <H1> or <H2> to verify that the matrices generated are not numeiic, but symbolic (algebraic) entities. 
Compute the algebraic outcome of the sum-product P by typing <P=H1 ' *H1+H2 . ' *con j (H2 ) ; >. Verify the algebraic values of the minors 
of P corresponding to P and P^ by typing P ( 1 : K/2 , K/2 + 1 : end) and P (K/2 + 1 : end, 1 : K/2 ) , respectively. Notice also that there are 
additional pairs of orthogonal vectors in the GABBA encoded matiices. This property will be exploited in section lTvl to devise a fully orthogonal 
decoder for the GABBA codes. In order to inspect the ^Masi-orthogonality property over real matrices, compute the the sum-product P with 
<P=H1 . ' *H1+H2 . ' *H2 ; >. Notice the interesting difference in the results. For example, for K = 8, the sum-product P is in fact diagonal! 
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Due to the block structure of Hx, it is clear that the noise independent- terms in equation ( I40t can be spht into 
the following independent components: 



K 



fK, = [[H«j"^T 0] 



f;,,= [[H«J°4 0] 



[0 [ulf-f] 



[0 [H^J°-]- 



(41) 



(42) 



It is left to prove that the matrices H^-j and which are referred to as equivalent encoded channel matrices, 
are identical. From equations (I24> and ( I26> : and from equations (125 > and (I26> . we have 



Ai 


Bi 


Bi 


-Ai 


B2 


-A2 


A2 


B2 



(43) 



(44) 



Substituting these into equations ( I41> and ( I42> . we have, straightforwardly: 



if 



-K2 



[AH B«]- 



[BH -AH]- 



Ai Bi 
Bi Ai 

Ai Bi 
Bi Ai 



B* 


-AS 




A* 


BS 





|Ai|2+|Bi|2+bt.b;+at.a;4h, 



(45) 



K 



BS 


-AS 




AS 


BS 





= |Ai|2+|Bi|2+Bi.BS+AT.AS ^ Hk/2, (46) 



which concludes the proof. 

Referring to equations J41t through (I44t . it is clear that equations ( 145 l l and ( I46l l can be combined into 

Hifi = H/f2 = H/f/2 = - (H"j - Ha'i + - H^J ' 



(47) 



which removes the need for partitioning the minors and H/f^ when computing the GABBA reduced encoded 
channel matrix. 

The reader is invited to inspect the ^Mflii-orthogonal decodabiUty of GABBA codes using the Matlab programs 
provided in Appendices HI and HIF . 

*To work with the examples of Appendices and |13 type <K=32;>. Generate the GABBA encoded channel matrix 
with the command < [ H, HI , H2 ] =GABBAEncodedChannelMat rix ( [ ] , K) ; >. Next, type <for k=l : K/2 , eval ( [ ' syms s' 
num2str (k) ] ) , end; > to define K/2 symbolic variables. Construct the symbol vector with its last K/2 entries reduced to zeros 
with <s=[];for k=l:K/2, eval ( [ ' s = [ s ; s ' num2 str (k) '];'])» end> and <szeros=[s; zeros (K/2, 1) ] ;>. Apply the 
GABBA encoder with <C=GABBAEncoder (szeros) ; >. Type <r=H* [szeros; conj (szeros) ] ; > to obtain the receive vector and, 
finally, type <rhat2=H ( : , K/2 + 1 : K) ' *r + ( r ' *H ( : , 3*K/2 + l : 2 *K) ) . ' ; > and <expand ( rhat2 ) >, obtaining a vector of algebraic 
zeros which confirms that the outcome of the linear combination given by equation 1421 is independent of the first K/ 2 symbols in sj^ ■ Similarly, 
type <szeros=[zeros (K/2, 1) ;s] ;> and <C=GABBAEncoder (szeros ) ; > to generate and encode a symbol vector with its first K/2 
entries padded with zeros. Compute the received signal vector with <r=H* [szeros; conj (szeros) ] ; > and the linear combination given 
by equation J^l with <rhatl=H ( : , 1 : K/2 ) ' *r + ( r ' *H ( : , K+1 : 3*K/2 ) ) . ' ; >. Finally, inspect the result with <expand (rhatl ) >. 
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In light of the similarities between equations ( I38> and ( I39l l and the equivalent system description given in equation 
( I16> . the matrix Ha/2 will be hereafter referred to as the GABBA reduced encoded channel matrix. 

Theorem |2l establishes that the GABBA codes admit a ^MOiZ-orthogonal decoder similar to that described in 
[27]-[29] and, therefore, can be seen as a complete generalization of the ABBA scheme previously known for 
limited number of antennas only. In the next section it will be shown, however, that a fully orthogonal decoder for 
the proposed GABBA codes also exists. 

IV. Orthogonal Decoding of GABBA Codes 
A. Properties of GABBA Reduced Encoded Channel Matrices 

The essence behind Theorem|2]is that the n-th order manifold structure of the GABBA encoded channel matrices, 
which ensures their ^'waiZ-orthogonality, ultimately enables the decomposition of the received vector into two vectors 
each dependent of mutually exclusively _ft'/2-tuples of transmitted symbols. Similarly, it will be shown that the 
symbol-by-symbol decodability of GABBA codes derives from additional orthogonalities in the reduced encoded 
channel matrix, beyond those proved in Theorem |2] 

To this end, let us first introduce the following result. 

Lemma 4 (Structure of GABBA Reduced Encoded Channel Matrices): 

GABBA reduced encoded channel matrices Ha-/2 in — l)-th order manifold matrices with generating function 
as given in equation M5\ . 

Proof: In similarity to Lemma|3] let us write the minors A^ and in equations ( I45> and ( I46t one step more 
explicitly, yielding Ai =i?(Ai,i, Bi,i), Bi = i/(Ai,2, 61,2), A2 = i/(A24, B24) and B2 = H(A2,2, B2.2). Then, 
the equivalent encoded channel matrix HAy2 can be written as 

|Ai,,P+|Bi,.|2+bt,.B^,,+AT_,;.a;,, AH,.Bi,.-B»,.Ai,,-BT^.A*2,,+ AT,.K2 
BH,.Ai,,-A» .Bi,,; + BT -A* -Al,.B* |Ai,,|2+|Bi,,|2+BT -B^.^+A^ -A*., 



A B 
B A 

(48) 

It is evident from equation J48t that Ha/2 has the structure described by equation (I25t . Proceeding in a similar 
manner it is straightforward to show that, furthermore, the minors A and B also have the same structure. In other 
words, Ha/2 can be constructed by iterating the generating function given by equation J25> n—1 times and, therefore, 
is a (71— l)-th order manifold matrix. ■ 

The fact that GABBA reduced encoded channel matrices Ha/2 have a manifold structure similar to that Ha" 
prompted us to investigate its orthogonality properties. One difficulty in dealing with the reduced matrices in 
algebraic form, however, is that its elements are sums of products of the channel estimates. A solution to this 
problem is to relabel the entries of Ha/2- A Matlab program that can be used to systematically relabel the entries 
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of GABBA reduced channel matrices is provided in Appendix IIII-AI *. Another solution is to simply rely on the 
proof of Lemma |4] and construct the relabeled GABBA reduced encoded channel as the manifold matrix of a 
channel vector of length with the . The corresponding Matlab code is provided in Appendix IIII-Bf °. Finally, one 
could also use the program given in Appendix lII-BI and generate the (n— l)-th manifold matrix as the second minor 
of the GABBA encoded channel matrix (i.e., the matrix <H2>), using a modified channel vector as input". 

Empowered by the ability to algebraically generate GABBA reduced channel matrices of any size, one can easy 
inspect their properties. Such an inspection quickly reveals that GABBA reduced channel matrices are real quasi- 
orthogonal (see definitions |6] and Q, i.e., there is a permutation that can be applied to a reduced GABBA encoded 
matrices such that the first half columns(rows) of the resulting permuted matrix are real-orthogonal to the remaining 
half columns(rows). 

For example, it is found that first and fourth columns(rows) of the 4-by-4 GABBA reduced encoded channel matrix 
are real-orthogonal to its second and third columns(rows). Likewise, the {1, 4, 6, 7}-th columns(rows) of the 8-by-8 
GABBA reduced encoded channel matrix are real-orthogonal to its {2, 3, 5, 8}-th columns(rows).. Analogously, the 
{1, 4, 6, 7, 10, 11, 13, 16}-th and {2, 3, 5, 8, 9, 12, 14, 15}-th columns and rows of the 16-by-16 GABBA reduced 
encoded channel matrix are mutually real-orthogonal, and so fourth. Though rather hard to prove in a concise and 
generalized fashion, these results are easy to verify algebraically (for example, using the Matlab programs given 
in Appendix IIIH . and indicate a clear pattern of real-orthogonality between the vectors of the GABBA reduced 
encoded channel matrices, which is conjectured to behave as follows. 

Conjecture 1 (Real Quasi-Orthogonality Pattern of GABBA Reduced Encoded Channel Matrices): 

Consider the function 

/ Llog2^'J-l \ 

Pn{x) = modf [\og^N\ - ^ + £ ' 2 J ' (49) 

where [xj denotes the largest integer no greater than x and mod(a;, 2) denotes the rest of division by 2. 
Let the sets Pg JV ™d Pi jv be respectively given by 

Po:N = {I7 • ■ • such that PN{i) ^Q,{i ^ N), (50) 

p^.j^ = {2, • • • , j, • • • }, such that pnU) = 1, (j ^ N). (51) 

'Before running the program, type in a value for the variable K, e.g., <K=32;> and generate the corresponding GABBA en- 
coded channel matrix and its minors with the command < [H, HI , H2 ] =GABBAEncodedChannelMatrix ( [ ] , K) ; >. Finally, type 
< [Hhat 1 , Hhat2 ] =RelabelH (HI , H2 , 1 ) > to generate and display the GABBA reduced encoded channel matrices. Notice that the 
outcomes <Hhatl> and <Hhat2> are generated independently and yield the same results, as proven in Theorem|2] It is also evident from the 
results that GABBA reduced encoded channel matrices are (n — l)-th order manifold matrices with generating function H{A, B) = , 
as proved in Lemma l4l 

'"To use the program, assuming that the variable K has been assigned a power-of-two value, simply type 
<H=RelabeledGABBAReducedEncodedChannelMatr ix ( 1 6 ) >. Notice that, in order to remain coherent to the notation used 
throughout the paper, the relabeled matrix outputted is of size K/2. 

"To generate the modified channel vector, type <h=[];for k=l:K/2, eval ( [ ' h= [h; h' num2str(k) '];']) r end> and 
<h=[h(K/4 + l:K/2) ;h(l:K/4) ]> into the command hne, after giving the variable <K> an appropriate (power of 2) value. 
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If Hjf/2 denotes the Kl2-by-K/2 GABBA reduced encoded channel matrix, then the permuted matrix given 
below is real ^'waii-orthogonal 

nK/2 = [[liKnY"'-'"^ [nK,2Y"-'<'-] ■ (52) 

That is, the pp.^-th columns(rows) of the A^-by-A^ GABBA reduced encoded channel matrix are real-orthogonal 
to its pj^ jy-th columns(rows). □ 
A Matlab program to generate the permutation indexes 7V ™d Pi.pj is provided in Appendix IIVI 
Conjecture [2 is algebraically verifiable for matrices of very large sizes and, therefore, is nearly as strong as 
theorem for most practical purposes. Indeed, using the programs provided in Appendices |III] and IIVI it is found 
that conjecture [2 is hold true for matrices as large as Matlab can handle in symbolic form'^, which is beyond the 
dimension of virtually any foreseeable MIMO application! The reader is invited to test our claims'^. 

B. Orthogonal Decodability of the Original ABBA Code 

The real q'Mfls/-orthogonality of GABBA reduced encoded matrices implies that the vectors r^i and Yk^, of 
equations J41> and (I42> . can be further decomposed into a total of 4 vectors, each dependent of mutually exclusively 
iC/I-tuples of transmitted symbols as follows. 

First, notice that the real q'Mai/-orthogonality property of 'H.k/2 ensures that, 

Ui 



2 2 



7^2 



(53) 



Let rjfj, rK27 Sj^/2 ™d 5^/2 be the permuted versions of the vectors yk^, yk2, Sk/2 and 8/^/2, respectively. Then, 
equations ( I41> and (I42l l can be rewritten as (see also equations ( I45> and ( I46» . 

rxi = 'Hk/2 ■ Ski 2-, (54) 
tk^ = 'Hki2 ■ hi/ 2- (55) 
Substituting equation ( I53> into these equations we obtain, 

(56) 
(57) 
(58) 
(59) 

'^This is dependent on the amount of memory available in platfomi used. We have verified the conjecture to hold over GABBA encoded 
channel matrices of sizes up to X = 1024. 

After generating an algebraic and relabeled GABBA reduced encoded channel matrix, construct the vectors Pq at ™d Pi jv 
with the command < [pO , pi ] =PermutationIndexes (1 : K/2) > and type <H ( : , pO) . ' *H ( : , pi) >, <H ( : , pi) . ' *H ( : , pO) >, 
<H (pO, : ) *H (pi, : ) . ' > or <H (pi, : ) *H (pO, : ) . ' > to verify the real-orthogonalities of the column- and row-wise partitions of <H>. 





= ['Hk/2 


rKi 


U:K/4 


= 7^1 • 


[sK/2 


U:K/4 


--Hi • 




IP0:K/2\ 




= [T^K/2 


rK^ 


D:K/4 


= 'H2- 


[sK/2 


D:K/4 


^n2- 




[Pl:K/2' 




— i'HK/2 


rK2 


U;K/4 


= ni ■ 


[Sk/2 


U:K/4 


= ni ■ 


h<] 


[P0:K/2- 


rK2,2 


= ['Hk/2 


■ rK2 


jD:K/4 


= n2 


[«A72 


jD:K/4 


= n2 


[sk 


[Pl:K/2 
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The last terms in equations (I56> through ( 15 9> are obtained by recalHng that the product H^^2 ' is a sparse 
matrix with the elements of Tii distributed across its P0 A72 columns and rows, and with the elements of 7^2 
distributed across its columns and rows, such that the half-partitions of Po if/2 respectively, the indexes 

of the symbols in the permuted vectors and and, likewise, Pi.K/2 is the set of indexes of the 

symbols in [si<-/2]'^'^"^ and [s/f/2]'^'^'*, respectively. 

In particular, for the case of i^T = 4, where the GABBA code reduces to the a variation of the original ABBA 
code, equations ( I56> through ( I59l l ensure the full decoupling of all four symbols encoded. In that case, it is trivial 
to show that Tii and ?<2 reduce to identical scalar numbers, such that equations ( I56> through ( I59l l need only be 
normalized accordingly to yield soft estimates of the encoded symbols. 

For the convenience of the reader, explicit calculations for the orthogonal decoding of the 4-by-4 GABBA code 
are put together in a Matlab script file provided in Appendix I V-AI One needs only to study and execute the file to 
verify the veracity of our claims. 



C. Orthogonal Decodability of GABBA Codes 

It was shown, so far, that the ^Mfl.s/-orthogonality of GABBA encoded channel matrices and the real quasi- 
orthogonality of GABBA reduced encoded channel matrices enable the orthogonal decodability of the 4-by-4 
GABBA code, which is a variation of the original ABBA scheme proposed in [27]-[29]. 

In this subsection, it will be shown that the orthogonal decodability is a general property of the GABBA codes. 
To this end, it is sufficient to show that the real g'Masi-orthogonality of GABBA reduced encoded channel matrices 
propagates to the minors T-Li and 7^2, and so forth. 

Conjecture 2 (Real Quasi- Orthogonality of Second-Order GABBA Reduced Encoded Channel Matrices): 

Let Til and 7^2 denote the minors of the real-product of two K/2-by-K/2 permuted GABBA reduced encoded 
channel matrices, as given in equation ( 15 3> . Then, the second-order GABBA reduced encoded channel matrix, 
defined as 

H^^4 = T-lI ■ 7i2 — 7^2 ' (60) 

is a real q'Mas;-orthogonal matrix. 

denote LIIC pClliiULCU VCISIUII Ul ±±^y4 



Furthermore, let the 7t|-/4 denote the permuted version of H?,, . such that 



4 4 



Such a permutated matrix is given by 



nl 

7^2 



(61) 



7i^/4 = [[n^^.r--^- [Hl./2]^--/-] , (62) 
where Pq.k/4 and ^1/^/4 are as in equations ( I50t and i5H . □ 
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Again, the veracity of conjecture |2] can be easily verified for values of K as large as Matlab is able to handle'"* 
and, therefore, has the strength of a theorem form most practical purposes. 

The similarity between Conjectures and 121 is evident. In fact, these results can be can be combined and 
generalized as follows. 

Conjecture 3 (Real Quasi-Orthogonality of n-th Order GABBA Reduced Encoded Channel Matrices): 
Let the n-th GABBA reduced encoded channel matrix defined''' by 

1 



H^V2 = 2 (HH^ • Hk, + . J , (n = 1) (63) 

H]j/2" = o<'r'y ■ nr' - mry ■ nr' (i< « ^ log.K). (64) 

Then, for all 1 ^ n ^ log2iir, the permutation 

^^72-^ = [[H^/2„]f--/- [H]J/2„]Pi-/-] , (65) 

is such that 



T ^T/n 
riK ■ rlK 





m 



(66) 



2" 2" 

with ?<" and l-Q reducing to a complex scalar at n = logjif. 

In other words, GABBA reduced encoded channel matrices of n-th order are real ^'wasj-orthogonal. □ 

The claims stated in conjecture |3] are, once more, verifiable to matrices of very large sizes and, therefore, are a 
solid results that can be used with the strength of a theorem for all cases of practical interest. A Matlab script file 
is provided in Appendix I VI- Al to facilitate this verification. 

The orthogonal decodability of GABBA codes is an immediate and obvious consequence of the real quasi- 
orthogonality of GABBA reduced encoded channel matrices. The following algorithm to obtain orthogonal estimates 
of all symbols encoded in GABBA encoding matrix can be inferred directly from conjecture |3l 

Algorithm 1 ( Orthogonal Decoding of GABBA Codes with a Multi-antenna Receiver): 

Consider a if-by-n,. matrix constructed from the collection of independent receive vectors YKy.ir ■ ■ tYk^-u^, 
each obtained at a different receive antenna, corresponding to the transmission of a A'-by-nt GABBA encoding 
matrix through a MIMO channel consisting of n^ transmit and n^ receive antennas. Assume that the n^-by-l channel 
vectors h„j.i, • • • ,h„j.„r at all receive antennas are perfectly known. Then, orthogonal estimates of all symbols 
encoded in GABBA encoding matrix are obtained as follows: 

'''For instance, input a value for the variable K, e.g., <K=32;>, and type in the command 
<H=RelabeledGABBAReducedEncodedChannelMatrix (K) ; > to generate a relabeled, algebraic GABBA reduced encoded 
channel matrix of desired size. Next, compute the second-order GABBA reduced encoded channel matrix algebraically through 
the real-product <HH=H.'*H;>, as in equation 1601 . Verify the commutativeness of the product as implied in that equation 
with the command <HH==H*H.'>, which should return a matrix of logical "ones". Then, compute the permutation indexes with 
< [pi, p2 ] =PermutationIndexes (1 :K/2) ; > and, finally, type in <HH(pl,p2)> and <HH(p2,pl)> to verify that the off- 
diagonal minors of the permuted second-order GABBA reduced encoded channel matrix are indeed null-matrices. If desired, type in 
<HH (pi, pi) > and <HH (p2,p2) > to inspect the diagonal minors corresponding to and T-tj in equation <611 . 

'^For convenience, equation 1471 is repeated in equation 1631 . 
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1) For every i-th receive antenna, construct the corresponding K-hy-2K GABBA encoded channel matrix as 
described in Lemma |2j 

2) For every i-th receive antenna, combine the received vector as shown in equations J41t and J42t . obtaining two 
independent reduced receive vectors tki A and rK2:i\ 

3) For every z-th receive antenna, use equation i63\ to construct the 1*^' order GABBA reduced encoded channel 
matrix ViK/2:i', 

4) Add the first reduced receive vectors Yki-a across the receive antennas, obtaining Yki:y.\ 

5) Add the second reduced receive vector Yk^a across the receive antennas, obtaining fA'2:s; 

6) Add the 1*" order GABBA reduced encoded channel matrix across the receive antennas, obtaining 11^/2" :s; 

7) Set n = 1; 

8) Multiply the aggregate reduced receive vectors yk^.t., (j = 1j • • ■ ■ 2") by the transpose of the aggregate n-th 
order GABBA reduced encoded channel matrix obtaining 2"+^ new aggregate reduced receive vectors, 
each dependent on mutually exclusive iir/2("+^^-tuples of transmitted symbols (see equations i56\ through i59l): 

9) Compute the permutation indexes Po:if/2" and Pi;k/2"'^ 

10) Obtain a higher-order GABBA reduced encoded channel matrices by partitioning 'Rk/2^:T,^ as in equation \65\ : 

11) If n < \0g2K — 1, increment n by 1 and return to step 8; 

12) Otherwise, the symbols have been fully decoupled. Normalize the combined signals'^ YKj-^, (j = 1, - • • ,2") 
by the GABBA reduced encoded channel matrix (which at this point has been reduced to a scalar) to obtain the 
soft estimate of the corresponding symbol. □ 

Notice that right after step 6 of algorithm ^ we have; 

nr nr 

YK^-.T. = E {yk^-.i + n^nO E (Hif/2:i • Sk/2 + ^^Kl■.^) = 'Rk/2:1: ' 8^/2 + n/f^^s, (67) 

1=1 1=1 

nr nr 

YK2:T. = i^K2:t + ^K2:i) = {^K/2:i ' 8^/2 + ^K2:i) = ^K/2:T. ' Sk/2 + ^Ki.T., (68) 

1=1 i=l 

where UK^-i, n^aii, yiki-.t, and nK2:T, are noise vectors. 

It is clear from equation (I48> that the diagonal elements of each I'*' order GABBA reduced encoded channel 
matrix Hx/2:i are real scalars given by X]j=i l^^j:! while all remaining elements are independent complex numbers. 
Consequently, the summation of the vectors yk^.t. and yk2:Y, across the receive antennas - which incurs in no loss 
of information, nor in any decrease of signal to noise ratio - significantly reduces the complexity the orthogonal 
soft-decoding procedure in a multi-antenna receiver, by eliminating the need for the parallel processing of the signal 
streams at all receive antennas independently. 

Notice also that the soft estimates obtained at step 12 of [2 are not in a trivial ascending ordered. In fact, 
although the symbols are split into two _fs72-tuples with indexes arranged in ascending order up to the first linear 
combination, at the following steps each of the tuples are split into two tuples arranged according to the sequences 

""The notation is slightly abused here for the sake of simplicity. 
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and Therefore, the soft symbol estimates obtained by the GABBA orthogonal decoder are ordered 

in accordance with the sequence computed as follows. 

Algorithm 2 (Order of GABBA Orthogonal Soft Estimates): 
Start with the matrix 



Pi = 



1 2 ■■■ K/2 

if/2 + 1 X/2 + 2 ••• K 



(69) 



Next, for n = {2, • • • ,log2-ftr}, partition the matrix P„ in half (columnwise) and stack the two partitions, 
obtaining 

[P ]P0:K/2" 

(70) 



where w ^^'^ Pi-n computed as described in conjecture Q 

Then, the j-th entry of the resulting column vector Piog^K is the index of the symbol corresponding to the j-th 
orthogonal soft estimate output by algorithm ^ □ 

A Matlab implementation of algorithms ^ and |2l combined into a fully orthogonal decoder for the GABBA 
codes, is provided in Appendix I VIII The reader is encouraged to test the generality, speed and numeric efficient of 
the technique'^. Try large values of K and over numeric symbol taken for different constellations to experience 
the remarkable low-complexity and flexibility of the GABBA decoder'*. Emulate the case of nj < K by simply 
puncturing arbitrary columns of the GABBA mother encoding matrix and padding the rows of the corresponding 
encoded channel matrix with zeros. 

As a final remark for this section, we emphasize that the decoding procedure described in algorithm \l\ is not 
to be confused with a zero forcing or any iterative technique such as the minimum mean square error (MMSE) 
decoder. Indeed, in the proposed decoder, all transmitted symbols are decoupled through linear combinations of the 
received signal vector, without any assumptions on the modulation scheme or any knowledge/estimation of the noise 
statistics, and without resorting to matrix inversions. Instead, as shown in the preceding sections, the orthogonality 
of the soft symbol estimated is ensured by the ^Mosi-orthogonality of the encoded channel matrix and the real 
g'Masi-orthogonality of the reduced encoded channel matrix of different orders. 

In fact, it is evident that the nested calculation described in algorithm can be reduced to a single-step 
implementation, where the linear combination weights Ffe:i and Fk:2 required to obtain the orthogonal estimate of 
each symbol Sk are first computed from the corresponding products of GABBA reduced encoded channel matrices, 

'^For instance, after giving tlie variables K and nt desired values, type <s= (randn (K, 1) +j*randn (K, 1) ) ; > to generate a symbol 
vector taking from a complex Gaussian constellation and and GABBA-encode it with <C=GABBAEncoder (s) ; >. Then, generate a Rayleigh 
block-fading channel matrix with the command <hh= (randn (K, nr) +j*randn (K, nr) ) ; > and compute the corresponding receive signal 
matrix with <rr=C*hh; >. Employ the GABBA decoder to obtain the soft orthogonal estimates of the transmitted symbols calling the function 
<sx=GABBADecoder (rr, hh) ; >. Finally, verify the results by listing the original and the orthogonal estimate symbol vector next to one 
another with < [s sx] >, or compute the norm of their difference with <norm(s-sx) >. 

'^Note that Matlab's floating point relative accuracy is of approximately 2-10^^''. Extremely large values of K may cause Matlab to overflow. 
This problem can be avoided by normalizing the intermediate quantities. 
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and then applied directly over the original received signal vector, as shown in equation ( 17 It 

Sfc F^^^-tk + r*j^-Fl2 = a/c • Sfc +nfc, (71) 

("r,nt) 

where ak is a real scalar proportional to ^ while fifc is the linear combination of the original noise 

terms rifc (see equation il6V ) with independent, complex coefficients. 

In other words, it can be said that the GABBA orthogonal decoding "algorithm" is a true extension of the 
Alamouti decoder, currently available for the case of nt — 2. The flexibility in handling with GABBA codes of 
different sizes (scalability), and the relatively easy expressions involved are, nevertheless, advantages of the nested 
implementation here described. 

V. Performance of GABBA Codes 

In this section, performance of GABBA codes in uncorrected fading channels of generalized statistics is studied 
both analytically and through simulations. 

Given the difficulty in computing explicitly the coefficients a^, as well as the corresponding coefficients multiply- 
ing the noise terms in rifc, the exact BER probability of orthogonally decoded GABBA codes is hard to determine 
in a generahzed and fashion. Therefore, we shall take the approach of deriving analytical expressions (bounds) 
corresponding to the ML decoder (whose complexity is prohibitive for large systems and high-order constellations), 
and compare the curves with simulation results obtained using the proposed GABBA orthogonal decoder. 

A. Performance of GABBA Codes with ML Decoding: Analysis of Exact BER Probability 

The BER probability of digitally modulated signals in a fading channel can be computed by averaging the 
corresponding error probability in the additive white Gaussian noise (AWGN) channel over the statistics of the 
fading process [35]. In a setting with nt transmit and rij. receive antennas, a ML-decoded rate-one STBC with full 
diversity is equivalent to a {nt x ?ir)-branch maximum ratio combined receive diversity system with receive power 
equal to n^, i.e., where the unitary transmit power is shared amongst all the nt transmit channels [36]. 

Consider an nt-hy-Ur MIMO set-up in which the constellation-energy-to-noise-power ratio (Es/No) of the ?i-th 
pair of transmit-to-receive antennas is denoted 7„, and where the envelope fading process corresponding to that 
diversity branch is described by a continuous probability density function (PDF) Pr„(7n)- Let Pb(7„) denote the 
BER probability of the digital modulation system over the n-th branch under AWGN conditions. Following the 
arguments in the previous paragraph, the BER probability of the full-rate full-diversity GABBA codes over this 
MIMO structure is given by 




(72) 
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1 ) M-ary PSK Modulation: 

The BER of M -ary PSK modulated signals over the AWGN channel were (believed to be) calculated exactly by 
Lee in [37]. A minor error (occurring at higher-order constellations) was recently discovered by Lassing et al. and 
corrected in [38], [39]. In particular, from [38, eqs. (2) and (3)], we have 

1 



A/-1 



Pb:PSKi-/\M) 



logjA/ 



J2 ^k - 



PSK • Pk:PS,Y.{n\M), 



k=l 



with 
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M 
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i=2 



(73) 



(74) 



where [x] rounds x to its nearest integer, and Pk:P5K{j\M) is the probability that a received M-ary PSK symbol 
falls k sectors away from the sector it belongs to due to the effect of AWGN, under the signal-to-noise ratio 7. 
A convenient formula for Pk-.psKillM) can be derived'^ from [40, eq. (24)], yielding 



Pfe:PSK(7|M) = 

ZTT 



e 




He) dO 



where, for uniform M-ary PSK constellations, we have 



5k- 



2k- 1 

M 
2k + 1 

M 



(75) 

(76) 
(77) 



Equation i75\ also appears (with a minor error, however) in [35, pp. 201, eq. (129)]. 

Substituting equation il5\ into equation ( I73> , and the consequent result into equation ( I72L one realizes that 
the moment generating function (MGF) approach (see [35]) can be easily applied in the computation of the exact 
average BER probability of M-ary PSK over fading channels. In fact, equation i72\ can be rewritten as 

1 



P 



b:PSK 



(7|M) 



21og2M 



fc=i 



with the functions gp^t^{S) and I{S,g,N) respectively given by 



HS,g,N) 



2Tr(l-(S) AT 
n=l 



9 



(78) 

(79) 
(80) 



TT b n=i ""V sin^(6')y 
where A^7„(a;) is the MGF of the PDF p^ijn)- 

The exact BER integral formula given in equation (I78> can be easily extended to the case of a STBC with rate p 
and transmit diversity order^** r/. To this end, one needs only to realize that in a linear STBC of rate p and diversity 
rj, each symbol is transmitted 1/ p times by each transmit antenna (so that the energy of the transmit constellation 



"To this end. notice that in our context [40, A* = 0]. From [40, eq. (18b)], this gives [40, A = -$]. Finally, let [40, A] be given by the 
left and right angles defining each PSK sector The difference between the formulas at adjacent delimiting angles gives equation <75l . 

^"Transmit diversity order is understood as the ratio between the diversity gain achieved by the STBC, to the number of transmit antennas. 
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increases proportionally); and that the diversity order rj < nt can be translated as a uniform loss in the diversity 
contributions of all branches in the systems. These observations lead to 



M-l 



^2 k=l 



\ ),nt ■ nr-lp,!^) ~ /('5j,.gpsK('5fc ),nt • n^lp,??) , (81) 



where 
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/((5,5,iV|p,??) 
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de. 



(82) 



Equation (18 U is a general, simple and exact formula for the BER probability of linear STBCs of rate p and 
diversity order 77 employing uniform Af-ary PSK modulation over an ensemble of uncorrected fading channels 
with arbitrary and unequal statistics. 

2) M-ary QAM Modulation: 

In similarity to the previous subsection, an exact formula for the performance of linear STBCs over Af-ary QAM 
constellation can be derived as follows. The exact BER probability for QAM modulation in the AWGN channel 
was derived by Cho et al. [41]. In particular, for regular Af-ary QAM constellations we have, from [41, eq. (16)], 

log2\/M 



Pi 



b:QAM 



(7|Af) 



1 



Pk:QAM{l\M). 



(83) 



In equation (I83> . -Pfc:QAM(7|Af) is the probability that the fc-th bit in the AZ-ary symbol is in error, which can 
be put from [41, eq. (14)] into the more convenient form shown below 
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where 



d.:QAM= (-l)L ^ 

and Q{x) is the Gaussian Q-function, given by 
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The similarities between equations J73> and ( I83> , and between equations ( I75> and (|84} are evident. Substituting 
equations ( I83t through (I86> into equation ( I72> . and invoking the MGF method, we obtain 
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Again, equation (I87> can be generalized to STBCs of rate p and diversity order 77 by replacing equation (|80| 
with equation (I82> . which gives 
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For convenience, a Matlab implementation of the above-derived analytical BER probability formulas over un- 
correlated and unequal fading channels of Hoyt, Rice, Rayleigh and Nakagami statistics [42]-[45] is provided in 
Appendix 

Closed form formulas corresponding to equations (18 U and ( I89> can most likely be derived (at least under certain 
conditions) using ideas found, for instance, in [35]. Given the simplicity and conciseness of the formulas provided, 
as well as the accurate results obtained with these equations, however, it is unlikely that the resulting expressions 
(which oftentimes involve terms dependent on factorials) would be of use. 

B. Performance of Orthogonally Decoded GABBA Codes: Simulations and Comparisons 

In this section, the performance of orthogonally decoded GABBA codes are compared to the exact formulas 
(bounds) derived above. 

1 ) Equipower Rayleigh Block-fading Channels: 

First, in figure [J the BER performance of a large GABBA code (rif = 64) simulated over independent and 
identically distributed (i.i.d.) Rayleigh fading channels is compared against the corresponding analytical BER 
probability bound. The simulated results (white markers) were obtained with the orthogonal, symbol-by-symbol 
decoder described in subsection IIV-CI while the analytical results (black markers) were computed using equation 
d87t . with p = 1 and ?7 = 1. 

A number of important conclusions can be inferred from the comparison. Fkstly, the remarkably low complexity of 
the scheme is emphatically illustrated. Indeed, in order to decode a single block, the ML-decoder, whose performance 
is associated to the analytical (analytical) curves shown in the figure, would require the minimization of a decision 
metric over all possible 64-tuples composed of 64-QAM constellation symbols, resulting in the prohibitive decoding 
complexity of 0(64^"'). Likewise, a ^MasZ-orthogonal decoding procedure would require the minimization of two 
independent decision metrics over all possible 32-tuples of 64-QAM symbols, yielding an equally prohibitive 
decoding complexity of 0(2 • 64^^). In contrast, the proposed orthogonal decoder requires merely the comparison 
of each orthogonal symbol estimate against the 64 possible symbols in the constellation, the same complexity of 
an uncoded system. 

Secondly, hard evidence is provided to the fact that the GABBA orthogonal decoder is not a mere zero-forcing 
technique. This is true because the performance loss of the symbol-by-symbol decoding procedure relative to its 
bound decreases with the number of receive antennas. In other words, no systematic penalty in signal-to-noise ratio 
is payed by the use of the GABBA orthogonal decoder It is difficult to determine exactly the reason why the 
simulated results obtained with the orthogonal decoder do not match exactly the bounds. The difference between 
the analytical and simulated curves shown in figure Q could, in principle, be caused by coloration of noise samples 
in the decoupled symbol estimates. Since the noise terms at the output of the decoder are linear combinations of 
the original noise vector, with coefficients given by sum-products of the uncorrected channel, noise coloration is 
rather unlikely, especially with large codes and more so when multiple receive antennas are used. This author is 
led to conclude therefore, that the problem is (even if partly, but substantially) caused by rounding errors inserted 
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by limitations of the simulation platform utilized (in this case Matlab). Anyhow, given the extreme contrast in the 
complexities involved by each of the alternatives compared, the results shown in figure^are nothing but remarkable. 

Finally, it is clearly shown that the diversity order attained by the GABBA STBC with its orthogonal decoder is 
the same as the one theoretically achievable with the ML decoder. This claim is consistent with the discussion of 
section nVl and is further supported by the results shown in figure |2] In this figure, the simulated and analytical BER 
probabilities of GABBA codes of various sizes over QPSK modulated signals in the Rayleigh fading channel are 
compared. It is found that even with a small number of receive antennas (n^ — 2), the performance of orthogonally 
decoded GABBA codes is close to the theoretical bound, and that the difference between does not increase with 
the number of transmit antennas fit, i.e., as larger codes are used^'. 

2) Unequal Block-fading Channels of Same Statistics: 

The point of figures ^ and 121 was to illustrate how the GABBA codes and the corresponding orthogonal decoder 
provide a realistic, high-performing solution for a large, flexible, linear STBC achieving fuH-rate and full-diversity 
over i.i.d. fading channels. 

Next, we turn our attention to the effects of unequal diversity branches and on the performance of GABBA codes 
or, in other words, to the effect of utilizing large GABBA codes in the presence of non-i.i.d. fading channels. Our 
objectives are two-fold, first, to verify that the impact of unequal power and/or unequal statistics of the space-time 
channels is less significant as the number of transmit antennas increases (larger codes); and, second, to iUustrate 
the flexibility and accuracy of the BER expression given in equation 

First, consider the case of channels with the same fading statistics but different powers [46], [47]. This case can 
be associated to a distributed scenario [12] with cooperative GABBA space-time encoded transmission, where the 
multiple cooperating transmitters are in the same environment (same channel statistics), but at different locations 
(different signal strengths at the receiver). The cooperation is assumed to be perfect, i.e., cooperating transmitters 
are perfectly synchronized and coordinated, such that their transmit signals arrive at the receiver at the same time, 
at every transmission epoch^^. 

It is also assumed that the average power of each transmit channel is an independent random uniformly distributed 
in the interval [0, Pmax]- Let the average power of the /c-th transmit channel be denoted by pk. We start by asking 
ourselves what is the power distribution profile (across the different transmit channels) that best describes the 
average scenario faced by the system. Notice that if the statistics of all channels are the same, all diversity branches 
can, without loss of generality, be relabeled into ascending order Mathematically we have: 

pi < • • • S% Pfc S% • • • 5% pK, (90) 

^'Except when going from nt = 2 to nt = 4. This is because at jit = 2 the GABBA code reduces to the Alamouti scheme, which admits 
a linear orthogonal decoder identical to the ML decoder. 

-^A practical example would be an indoor OFDM multi-cast system, where several base-stations are reliably connected to one another 
(e.g., via cable) and transmit the same signal (e.g. a video steam) cooperatively to a wireless user. The OFDM set-up not only provides for 
memoryless channels but also ensures that a synchronous transmission (synchronized acess-points) results in synchronous reception for all users 
in the premises. 
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The random pk is referred to as the fc-th order statistic of the ensemble {pi, ■ ■ ■ ,pk}- For P^ax — 1, the expected 

value of the fc-th order statistic pk is given by [48, pp. 665, eq. (7.9.20)] 

1 

P.= (^^iyy(^7-^--(i-)--rf- (91) 



The integral in equation i9H is a variation of the well-known Beta function [48, pp. 544, eq. (6.12.1)] (a.k.a. 
Euler's integral of the first kind), which has no general solution. Fortunately, for {k, K} E N, it is easy to show 
(see Appendix Hx} that 

.^■n-.r-^d.^^-^^. (92) 



From equations ( I9U and i92\ . we finally obtain, for a general value of P^^x, 

Pk = ]f^^max. (93) 

In figure|3j the simulated and analytical performances of GABBA codes of sizes 4 and 64, over 16QAM-modulated 
symbols in an ensemble of Rice block-fading channels with equal and unequal (linear) power distributions are 
compared. It can be seen that the impact of an unequal power distribution on the ergodic performance of GABBA- 
encoded systems loses relevance when large codes are used. As a side product of the comparison, it is once again 

verified that orthogonally decoded GABBA codes (simulated results) perform very close to the analytical bound 

2k 

(in this case given by equation (|^ in combination with equation J82> . with N — K, p — 1 and 7„ — —). 

K 

3) Block-fading Channels of Different Statistics: 

Next, the performances of GABBA codes in uncorrected block-fading channels with equal and unequal statistics 
are compared. Results obtained both analytically and through computer simulations are shown in figure|3 The curves 
corresponding to the unequal fading channels were obtained by assigning to the fc-th branch a fading severity factor 
ruk and a power factor ilk such that nik spans the interval [0.5,4] uniformly, with nik+i > nik, while ^17^ = 1, 
with ilk+i < ^k- In other words, the channels experiencing more severe fading are also the most powerful ones, 
such that no particular channel in the ensemble is "dominant". This distribution of parameters is selected so as to 
prevent a few channels to dominate the others, ensuring that, instead, all branches have significance in providing 
diversity. This selection of parameters is associated to scenario faced by peer-to-peer MIMO systems employing 
diversity-optimized multi-mode adaptive antennas with orthogonal pattern adaptability^^. 

The comparison in figure |3 demonstrates how the use of a large GABBA code results in more robustness against 
shadowing and other effects that may affect the power distribution and fading statistics of the diversity branches 
available in peer-to-peer MIMO systems with multi-mode antennas. 

^^It is assumed here that the beamforming strategy is such that narrow beams are pointed towards stronger multipath components while wider 
beams are formed towards portions of the space characterized by more dispersive multipath clusters. A consequence of such approach would 
be that the narrower beams would collect less energy (fewer multipath components), but exhibit less severe fading, In contrast, while the wider 
beams would collect more energy (more multipath components), but also exhibit severer fading. 
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The unitary rate and full diversity achieved by GABBA codes ensure that they outperform all (presently known) 
alternative STBCs of linear construction and comparable complexity (orthogonally decodability). It is, therefore, of 
little use to show comparisons between GABBA codes and other linear STBCs such as the Alamouti scheme [2], 
the OSTBC [3], [4], and the maximal-rate minimum-delay (MRMD) codes [22]. The comparison of these codes 
can, nevertheless, be easily performed analytically utilizing equations (I78> and ( I87> . by simply selecting for each 
technique the corresponding values of p (for rate) and r/ (for diversity order). 



In this section, the performance of the proposed GABBA codes is analyzed from a capacity point of view. 

The Shannon capacities (over Gaussian constellations) of MIMO systems in ergodic and deterministic fading 
channels were derived in [16], [30]. More recently, the modulation-constrained rate capacity (maximum average 
throughput) of MIMO systems over Af-ary PSK constellations has also been derived [18]. It is known, however, 
that the rate capacity of linear STBCs (which do not provide multiplexation gain) is bound by the Shannon capacity 
of the corresponding SIMO channel [17]. 

It is of interest, therefore, to compare the average throughput theoretically achievable by GABBA codes in 
systems utilizing different modulation schemes over ensembles of nt-hy-rir antennas against the ergodic capacity 
of the equivalent l-by-n^ SIMO channel. Given the exact BER probability formulae derived in section IV-AI a 
simple yet effective approach at hand is to consider the rate capacity of GABBA codes with hard-decision on the 
orthogonal symbol estimates. In this case, one can simply considering the binary symmetric channel (as seen from 
the input of the space-time block encoder up to the output of the corresponding decoder, with the physical MIMO 
fading channel included), with the transition probabilities given by Pi, and 1 — Pi,, respectively [49]. This leads to 



where is computed with equations ( 17 8> or ( I87K at the given signal-to-noise ratio, considering the corresponding 
number of transmit and receive antennas rit and n^, and with parameters p = I and rj = 1^^. 

Using equation ( I94> . the hard-decision modulation-constrained rate capacities of GABBA codes of different 
sizes, achieved with different modulation schemes, in both dominant transmit (n* > n^) and receive (rir > rit) 



found that GABBA-encoded systems with optimized modulation comes to roughly 1 bit/sec/Hz from the Shannon 
SISO capacity. This is a rather remarkable performance when considered that the Shannon capacity is related to 
an unfeasible Gaussian modulation and soft-decision decoding. In contrast, the modulation-dependent rates shown 
for the GABBA scheme can be reached at significantly lower complexities, through concatenation with efficient 
inner codes such LDPC or Turbo Codes [50], and employing hard decision at the symbol-level. In other words, it 
is learned from figures |5(a)| and |5(b)| that the considerably more complex soft-symbol decoding of GABBA codes 
(concatenated with an inner code) would only improve the rate of the system by 1 bit/sec/Hz. 

^''For simplicity, all these parameters are omitted from tlie notation of Pj,. 
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diversity scenarios, are compared against the SIMO Shannon capacity in figures |5(a)| and |5(b)| respectively. It is 
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It is also noticeable from figures |5(a)| and |5(b)| that for lower signal-to-noise ratios and relatively lower-order 
constellations (up to 64QAM), the GABBA codes considered {rit = 4 and rit — 16) actually come to less than 
1 bit/sec/Hz from the Shannon bound. This result is further illustrated in figure |6l where the achievable rates of 
balanced MIMO systems (nt — rir) employing optimum modulation and GABBA codes of different sizes^^ are 
compared to the corresponding SIMO Shannon capacities. It is interesting to notice that MIMO system with the 
GABBA scheme and optimum modulation exhibits the same general behavior as the "ideal" Shannon bound when 
it comes to rate increase in proportion to the number of antennas used in the system. 

VII. Conclusions 

In this paper, a family of linear space-time block codes (STBCs) of unitary rate and full diversity, systematic 
constructed over arbitrary constellations for any number of transmit antennas are introduced. The codes were 
obtained from a generaUzation of the existing ABBA STBCs, a.k.a ^Mosi-orthogonal STBCs (QO-STBCs) and, 
therefore, are denominated generaUzed ABBA (GABBA) codes (see section III-Bt . 

It was also shown that the GABBA codes admit fully orthogonal (symbol-by-symbol) decoding. Two key ideas 
underlie the derivation of the orthogonal GABBA decoder The first is a transformation applied to the system 
equation (see section El-AL that exposes the relationship between the code structure and the interference amongst 
symbols embedded in the same encoding block, characterizing it in terms of an equivalent encoded channel matrix 
(see section IIII-B> . The second is the exploitation of partition orthogonality properties present in structure of the 
GABBA codes, which are shown to be transferred to the corresponding encoded channel matrices. It was shown that, 
as a consequence of these techniques/properties, the received signal vector can decomposed into lower-dimension 
tuples, each dependent only on certain subsets of the transmitted symbols, via simple linear combinations with 
coefficients derived directly from the channel estimates. Orthogonal decodability results from the nested application 
of such linear combinations, with no matrix inversion or iterative signal processing required. 

The performance of the proposed GABBA codes were studied both analytically and through computer simulations. 
The exact BER probability of the codes with ML the decoder over generalized (uncorrected) fading channels was 
derived analytically, under the assumption of perfect channel knowledge at the receiver (see section IV-A> . The 
formulae were then compared to simulation results obtained with the proposed orthogonal decoder The comparison 
reveals that the proposed GABBA solution, despite its very low complexity, achieves nearly the same performance 
of the bound corresponding to the ML-decoded system, especially in systems with large numbers of antennas (see 
section |V-B> . 

A brief analysis of the average transmission rate achieved with GABBA codes was also provided. The analysis 
indicate that the transmission rates achievable with the full(unitary)-rate full-diversity GABBA codes over optimized 
constellations, with the corresponding orthogonal decoder and hard symbol decision, have the potential to come 1 
bit/sec/Hz from the theoretical bound given by the Shannol SIMO channel capacity. It has been shown that higher 

^^The curves are the envelopes of all possible regular modulations ranging from BPSK to 4096QAM. with hard-decision at the symbol level. 
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rates can be achieved using unitary-rate STBCs as building blocks in a layered structure. Such approaches have 
already been demonstrated for the known ABBA scheme and OSTBCs [51, Chap. 9], [52], [53]. It is foreseeable, 
therefore, that these techniques extend naturally to the family of GABBA codes here introduced, with the added 
flexibility, scalabiUty and low complexity to its advantage. 

In addition to the equations, proofs and figures found in the main body of the paper, extensive material is provided 
in the appendices to allow the verification of all claims, including expUcit examples and Matlab implementations 
of the encoder, the decoder and the BER formulae. 



February 1, 2008 



DRAFT 



29 



GABBA Codes in Rayleigh Fading 
Q (n^ = 64, 64QAM) 

1 t I ! ! I I j j ! 




SNR(indB) 

Fig. 1. Performance of the GABBA code C64 over 64-QAM modulation in Rayleigh block-fading channels with various numbers of receive 
antennas. Simulated results (white markers) were obtained using the orthogonal GABBA decoder described in section llV-CI while analytical 
results (black markers) were obtained with equation <87L with p = 1 and r} = 1 (full-rate full-diversity). 
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Fig. 2. Performance of GABBA codes of different sizes over QPSK symbols in Rayleigli block-fading channels with 2 receive antennas. 
Simulated results (white markers) were obtained using the orthogonal GABBA decoder described in section llV-CI while analytical results (black 
markers) were computed using equation 1871 . with p = 1 and rj = 1 (full-rate full-diversity). 
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GABBA Codes in Rice Fading 
(n =4, m = 1.5,J6QAM) 
1 ^ 
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^ Eq. (n^ = 64) 
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Fig. 3. Performance of GABBA codes of different sizes over 16QAM symbols in Rice block-fading channels with 4 receive antennas. Simulated 
results (white markers) were obtained using the orthogonal GABBA decoder described in section llV-CI while analytical results (black markers) 
were computed using equation 1871 . with p = 1 and r] = 1. Both simulated and analytical results are given for channels with equal and unequal 
(linear) power distribution profiles across the transmit channels (see section IV-B.2t . 
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GABBA Codes in Generalized Fading Channels 




10 11 12 13 14 15 16 17 18 19 20 10 11 12 13 14 15 16 17 18 19 20 

Es/No (in dB) Es/No (in dB) 



Fig. 4. Performance of GABBA codes of different sizes over 8PSK symbols in generalized block-fading channels with 2 receive antennas. 
Simulated results (white markers) were obtained using the orthogonal GABBA decoder described in section llV-CI while analytical results (black 
markers) were computed using equation <87l . with p = 1 and r] = 1. Both simulated and analytical results are given for channels with linear 
power distribution (see section IV-B. 21 and fading severity factors m spanning the interval [0.5,4], as well as equipower channels, and a mean 
fading severity factor m = 2.5. 
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(a) Dominant transmit diversity. 
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GABBA Throughput x Shannon Capacity 
(m = 1;nj = 4; n_.= 16) 
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(b) Dominant receive diversity. 



FiSrutirv jVS^W^'^'™ constrained achievable rate (ttiroughput) of GABBA codes of different sizes in Rayleigh fading channel and with dj|jggBf 
modulation schemes against signal-to-noise ratio. 
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GABBA Throughput X Shannon SIMO Capacity 
(m = 1 ; n^ = n^; SNR = {0.3.6.9} dB) 
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Fig. 6. Throughput of GABBA codes against number of transmit and receive antennas (nt = rir) in Rayleigh fading channel with optimized 
modulation. 
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Appendix I 
GABBA Encoder 



A. Example: GABBA Mother Encoding Matrix (K = 32 j 



C32 



siO sil ^12 S13 ^14 



si0^15^16 ^13 ^14^19- 



^21 



322 



-S9-«10-*'ll - 

Sll S12 -69- 
-**12 ^11 sio 

S13 ^14 S15 
-^14 ^13^16 
-^15^16 313 

3J6^J5-«14 

-^18^17 -'go 

~*j0"*j9~^j8" 
""91 "53- 
■^J2^J1 »J4 
^J3 =J4^21 
^J4-^J3-«22^ 
-^J5 =J6 =27 - 
"^J6^J5 °28 
-^J7 =J8^25 
-^28-^27^26- 

^g2~*'g9 
~**32~**31~*30" 



-^8 31 32 
-^7 ^2 31 
■^6 ^3 ^4 
-^5 34 ^3 
^12^13^14- 



«13 
"16 - 

315 

310- 

^9 



'go =gl-»g2 
'19 =g2 =gl 
'18 ^g4 
>J7 =g4 =g3 
'i4^J7 =18 
'83-18^17 
■82^J9 =20 
'gl^gO^lO 
'g8 =g9^J0 
'g7 =50 =29 
'g6 =51 ^§2 
'g5 =^2 =51 
'32 ^g5 =g6 

'ji^ge^gs 

'jO^gT 'gS 
'29~®28 ~*27 



^9 - 
310 
"g3 
"g4- 
3gl- 
3g2 
19 " 



34-«13- 
33 =14- 
32 315 
31 ^16 
"16 =1 
«15 -32 
«14 -33 
"13 "4 



=18 =31 
=17 =52 

= gl ^18 
= 50 ^19 
= g9 ^20 
= g8^gl 
= g7^g2 
"g6^g3 
*25-*24 



'13 =16 
'16 ^13 
'15 =14 



^4 

-=3 



"g6^g7 
= g5^g8- 
=g8 =25 
=g7 =g6 
"30-51 
= 29 —§2 
=g2 =29 " 

=31 =|o 

18 19 
=17 =gO 
= gO-J7 
= 19 —18 
=g2 =g3" 

^Si^S-* 

= g4-gl 
*23— 22" 



= 9 =16- 
= 16 =9 
= 15 —10 
= 14—11 - 
= 13 =12- 
=4 =5 
=3 —6 
=2 —7 



■g8-g9 

'g7-3o- 
'ge-ji 

■g5-32- 
'§2 =25- 
'31 =§6 
'30 =27 
'29 =28 
'20 =21" 
'19 =22 
'^8 =23" 
'17 =§4 
'g4-l7 
'g3-18- 
>g2-J9 
*21— 20- 



15—14 =13 =20 - 
10 
= 9 - 



'11 =12—21 
'12 =11 =22 —21 



= 9 
"10 

= 7 
—8 

"5 
—6 



10 =23 
= 9—24 
= 8—25 
= 7 =26 
= 6 =27 
= 5—28 
= 4 =29 
= 3 —30 =29 - 
= 2 —31 



23 
26 
—25 
"28 
= 27 



330 



= 21 —22 
322 —21 

— 27 —28 
= 28 —27 

— 25 —26 
"26 —25 



331 



332 



332 =31 



'52 =51 =g 

'g9 =J0 =3 

10 - go =J 

'g7 =g8 =5 

'g8-g7 =g 

'g5-g6 =7 

16 =g5 =5 

13 =g4 =9 

'g4-g3 "10 

11 - g2 =11 
12 =21 =12 
'j,9-g0 =13 

10 '19 

■J,7 "3,8 %5 

'18—17 =16 



= 3g 
-3 

'I 



=22 =21 
= 23 —24 
= 24 —23 
= 17 =18 
■=18 =17 
= 19 —20 
= 20 —19 
«29 —30 
"30 —29 
"31 =32 
«32 =31 
= 25 —26 ■ 
= 26 —25 
= 27 =28 - 
=2g =2J 



= 23 =24 =25 =26 =27 =28 =29 =30 =31 =32 
= 24 =23 —26 =25—28 =27—30 =29—32 =31 



= 21 =22 —27—28 =25 =26—31 —32 =29 



= 22 =21 



= 19 
= 20 



20 
19 
18 

— 18 =17 

— 31 —32 
"32 —31 

—29 —30 
"30 —29 



328 - 

329 - 

330 - 
331 
332 
317 
'18 
'19- 
320 - 



= 26 
= 31 - 
= 32 - 
= 29 - 
= 30 - 
= 19 
= 20 
= 17 
= 18 



= 27 —28 —21 —22 —23 - 
= 28 —27 =22 —21 =24 " 
= 26 =23 =24—21- 



""Ig "11 °12 " 
= 9 -Ig -Jl - 
-12 J9 -Ig - 
=11 =10 J9 - 

^16 -15 
""16 '13 -J,4 
"15 =14 =13 



2g -25 
= 1 -§ 

=8 =1 



i 
I 
=2 



■13 =14 
'14 — j,3 
'15 =16 -JS =14 
'Ig "15 -14 -13 
°.9 -IS -Jl =12 

'1° ^9 "IS "11 
'J,l -J2 %9 -Ig 
'12 =11 "10 =9 



'10 

'11 

'12" 

'13 

'14- 

'15 
'Ig- 



'25 =32 
'32 =25 =26 
'31 —26 =25- 
'30—27—28 
'29 =28-27- 
'20 =21 =22 
'19—22 "21- 
'18—23—24 
'17 "24-23- 
'24 "17 "18 
'23—18 "17- 
'22—19 —20 
11 =50-19" 

'19 n5"j6 

"9 'Ig '15 

'16 —9 =lg 

'15—^0 -^9 

;i4-Jl =12 

"5 4§ 
"s "J 3 



30 



330 =29 

=27 =28 

= 28 " 

=25 =26 

=26 =25 

=23 =24 

=24 "23 

=21 "22 

=22 "21 

=19 =20 



'18 =17 
'J5 '16 
'16"15 
'13-14 
'14 =13 



'12 
= 11 
= lg 
"S 



"J "8 



=5 



g 5 
"a '1 
"4 "a 
!i "s 



B. Matlab Function: Generator of GABBA Mother Encoding Matrices 
function C = GABBAEncoder (s) 

% Input: vector s of numeric or symbolic entries corresponding to symbols 
% Output: GABBA mother encoding matrix C 

% Ex.: syms si s2 s3 s4 s5 s6 s7 s8 s9 slO sll sl2 sl3 sl4 sl5 sl6 

% s = [si s2 s3 s4 s5 s6 s7 s8 s9 slO sll sl2 sl3 sl4 sl5 sl6] ; 

% C = GABBAEncoder (s) ; 

k = length (s); % Must be a power of 2 

s = reshape (s, [1 1 k] ) ; 

while k > 2; 

for n = 1 : k/2 , 

C ( : , : ,n) = ABBAl (s ( 

end 
s = C; 
clear C; 
k = k/2; 

end 

C = ABBA2 (s(:,:,l),s(:,:,2) 

% Auxiliary functions - 

function C = ABBAl (x,y); C 
function C = ABBA2(x,y); C 



2*n-l),s(:, :,2*n)); 



[ [x y] ; [-y x] ] 
[ [X y] ; [-y' x' 
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Appendix II 
GABBA Encoded Channel Matrix 



A. Example: Minors of GABBA Encoded Channel Matrix (K ^ 32) 



H 



32i- 



hi 


h2 


'•3 


h4 


h5 


hg 


h^ 


hg hg 


'•10 


hii 


'•12 


'•13 


hi4 


'•15 


'•16 


hiQ h20 h2i h22 '^23 ^24 '^25 ^^26 ^27 '^28 ^129 '^30 ^131 '^32 


'•2 


-hi 


''4 


-^3 


hg 


-h5 


'•8 


-h-jhia 


-hg 


'•12 


-hii 


hi4 


-'•13 


'•16 


-'•15 


'^18 -'^17 '^20-'il9 '1^22 -^121 '1^24 -'^23 '^26 "'^25 ^^28 "'^27 ^^30 "'^29 'i32-'i31 






-hi 


-h2 


'•7 


hg 


-hr, 


-hghii 


'•12 


-hg 


-'•10 


'•15 


'•16 


-'•13 


-hi4 


'^19 '^20-'^17-'^18 ^^23 ^^24 -''^21 -''^22 '^27 ^^28 -^^25 "'^26 ^^31 ''^32 "'^29 "'^30 


'•4 


-h3 


-h2 


hi 


'•8 


-h7 


-hg 


h5hi2 


-hii 


-'•10 


hr) 


'•16 


-^15 


-hi4 


'•13 


^120 -'^19 -'^IS '^17 '^24-'^23-''^22 2 1 '^^28 "'^27 "'^26 ^^25 ^^32 "'^31 -^^30 ^^29 


/i5 


hg 


h7 


ft' 8 


-hi 


-h2 


-h3 


-'•4 '•13 


hi4 


'•15 


'•16 


-hg -hio -hii -hi2 


^21 '1^22 ^^23 '^24 -'^17-'il8 -'^19 -'^20 '^29 ^^30 ^131 '^32 -^125 -^^26 -^^27 "'^28 


he 




''8 


^7 


-^2 


hi 


-h4 




-'•13 


'•16 


-^15 


-'•10 


hg 


-'•12 


hll 


'^22-'^21 ^^24 -'^23 -'^IS '^17 "'1^20 ^19 ^^30 "'^29 ^^32 "'^31 "^^26 ^^25 -^^28 ^^27 


h7 


''8 


-h5 




-h3 




hi 


h2hi5 


'•16 


-'•13 


-hi4 


-hii 


-'•12 


hg 


'•10 


^^23 '^24 -'^21 -'^22 -^^19 -'^20 '^17 ^118 '^31 ^^32 -^^29 "'^30 -^127 "'^28 ^^25 '^26 


hg 


-h^ 


-^6 


'•5 


-h4 


'•3 


'•2 


-hihig 


-^•15 


-hi4 


'•13 


-hi2 


hll 


'•10 


-hg 


^^24 -^^23 -'^22 '^21-'^20 '^lO ''^ 18 "'1^17 '^32 "'1^31 -^^30 '^29 "^^28 '^27 ^^26-^^25 


hcj 


'"10 


hii 


'•12 


'•13 


hi4 


'•15 


hig -hi 


-h2 


-h3 


-h4 


-h5 


-hg 


-h7 




^^25 '^26 ^^27 ^^28 '^29 ^^30 '^31 ^^32 -^117 "'^18 -^^19 "'^20 -^^21 "'^22 "'^23 -^124 


'"lO 




'•12 


-hii 


hi4 


-'•13 


'•16 




hi 




'•3 


-hg 


'•5 


-^8 




^26-^25 ^28^27 ^■?,Q^29 32 ^31 ^17^20 ^19^22 '^21 "'^24 '^23 


'"ll 


'"12 


-hQ -hio 


'•15 


'•16 


-'•13 


-hi4 -h3 




hi 


'•2 


-h7 


-h8 


hr, 


h5 


^^27 '1^28 -'^25 -^126 '^31 ^^32 "'^29 "'^30 "'^ 19 "'^20 ^^17 '^18 "^^23 "'^24 ^^21 '^22 


''12 




-''10 


hg 


'•16 


-^15 


-hi4 


'•13 -'•4 


'•3 


'•2 


-hi 


■*8 


hr 


hg 


-hr, 


^^28 -^^27 -'^26 ^^25 '^32 "^^31 "'^30 ^^29 -^^20 ^^19 ^^18 "'^17 -^124 '^23 ^^22 -^121 


'"IS 


'114 


'•15 


'•16 


-hg -hio -hii -fii2 -hs 


-hg 


-h7 


-^s 


hi 


h-l 


'•3 




^^29 '^30 ^^31 ^^32 -'^25 -'' 26 -'^27 -'^28 -^^21 -'^22 -^^23 -'^24 ^^17 '^18 ^^19 '^20 


'114 


-'"13 


'•16 


-hl5 


-'•10 


hg -hi2 


hii -hg 




-h8 


h7 


'•2 


-hi 




-^3 


^30^29 ^^32 -^^31 -'^26 ''^T) "'^28 ''■27 "^^22 ^^21 -^124 '^23 'il8-'^17 ^20^19 


''15 


'il6-'*13-^14-'ill-'il2 


hg 


'•10 -^7 


-'•s 


'•5 


hg 


'•3 




-hi 




''•32 -'^29 "''■30 "'^27 "''■28 ^^25 ''■26 "''■23 "^^24 ''■21 ^^22 ''■19 ^^20 "'^IT "''■18 


'*16-'*15-'*14 


'•13-'*12 




'•lO 


-hg -hs 


h7 


'•6 


-''5 




-h3 


-h2 


hi 


''■32 "'''31 "^^30 ''■29 "^^28 '''27 '^26 "^^25 "'''24 ''■23 ''■22 "^^21 '''20 "'^IQ "^^18 '^17 



H; 



322- 



17 18 -''19 ''20 -''21 
''18 '' 17 -''20 -'' 19 -''22 - 
''19-''2() ''IV -'' 18 -''23 



> ''26 ''■27 -^'28 
■,-^'2:, ''28 ''27 
7 ''28 -''25 ''26 



'20 



''19 



'■18 



' 17 



21 - 



'■'21 - 



'■28 



■ -/i 



''21 -''22 ^'23 ''24 ''IV 
''■22 ''■21 -''24 -'"'■23 '^18 
'^23^24 'i21-'^22 '^19 
'^24 '^23 '^22 '^21 '^20 
''•25 "'''26 "'^27 '^28 -'^29 
''26 ''•25 -''■28 -'^27 -'^30 
''27 "''•28 ''■25 -^26 -'^31 
''28 ''■27 ''■26 ''25 -'^32 
''29 -''■30 -''■31 '^32 '^25 
''30 ''29 -''■32 "''•31 ''■26 
''•31 "'^32 '^29 "''30 ''■27 
'^32 '^31 '^30 '^29 '^28 



20 -^'29 



L ^'32 

''^17 -''^20 -'^19 -'^30 -'^29 ''^32 ''31 
-/i20 '^17 -''■18 -''■31 '^32 -''29 '^30 

hiQ his '^17 -'^32 -^31 -'^30 -'^29 

'^30 '131 -'1 32 '^17-'^18-'''19 ''■20 
-''29 '^32 ''31 ''18 '^17 -'^20 -'^19 

''32 -''29 ''30 ''19 -''■20 ''17 -''18 
-'^31 -'^30 -''29 ''20 ''19 '^18 ''17 
-''26 '^27 ''28 ''21 -''22 -''23 ''■24 

''25 "'^28 -''■27 ''22 ''21 -''24 -''■23 
-''28 ''■25 -''■26 ''23 -''24 ''21 -''22 

''27 ''■26 ^^25 ''■24 ''23 ''22 ''21 



'32 



■31 



''29 



-''25 ''2(i ''27 -''28 
-''■26 -''■25 ''■28 ''■27 
-''27 ''■28 -''25 ''26 
-'^28 -'^27-''26 -'^25 
-''21 ''■22 ''23 -''■24 
-''22 -'^21 ''24 ''23 
-''23 ''■24 -''■21 ''22 
-'^24 -''23 -''22 -''21 
'^17 -'^18 -''19 '^20 
his ''17-^20-^19 
'^19-'^20 ''17"''18 
''20 ''19 '^18 '^17 



-''1 h2 

-412 ^'1 

^'3 ''4 

^'-1 ^^3 

^'5 ''6 

-h,5 -hs 

-h,7 hg 

-h,s -h7 

-hQ hiQ 

-hio -''9 
-hll ''12 
-''12 -''11 - 
-hi3 hi4 

-''■14 -'^13 
"'^15 ''■16^ 
-'^16 -'^15 - 



-III h2 

^'2 ^■l 

''V -*'« 

hg h-j 

-flQ hQ 

-hQ -h^ 

''■11 -'^12 

hi2 'ill 



''5 ^'6 
hg h5 
hy -h,g 



-''lO -''■11 

''9 -''■12 - 
-''12 ''■9 " 
hll ''■10 



hi 



-,-'M 



^' I 



/'9 



-4i2 -hi 
-/13 h4 
-h4 -h3 

''13 -''14 
hi4 hi3 



h4 h3/ti4 hi3-/ii(5 

"''1 ''2 ''15 
-h2 -hihie 



-!i< 



10 



-hi 



'10 -''9 
'15 "''16 
'16 '^15 
'13 '^14 



'14 - 



^13 



hie '^15 

-hg h-iQ 
-hio -^9 
-hll '^12 
-hi2 -/i-ii 



-''16 -''15 ->T-2 
''13 -''14 -''3 
''14 ''■13-^4 
''11 "''■12 "''5 
hi2 hii-hg 
-hg hiQ-hy 

"'^lO "''■9 "^8 



-^■16 "13 

h 1 5 h 1 4 

''2 ''3 

-hi h4 

h4 -?ii 

-h3 -h2 

hg hj 

-h5 hg 

hg -h5 

-hj -Hq 



'■16 

''16 ''9 

-hi5 hlQ 

-hi4 hll ■ 

''13 '^12 

h3 hQ 

-hi hg 

-hg -hi 

hy -h2 

hg -h3 

— hg — h^ 



11 "1 

-'•13 '•I 6 '^15 
'•16 -''13 ''14 
— h 1 tr^ — h 1 4 — h 1 

-''10-''ll ''12 
''9 -''12 -^11 



^12 



"10 



hll hio ''9 

-hg -h7 hg 

h 5 — h g -h 7 

-''8 ''5 -''6 

h7 hg hg 

''2 ''3 -''4 

-hi h4 hg 

h4 —hi h2 

-''3 -*^2 -^1 



B. Matlab Function: Generator of GABBA Encoded Channel Matrices 



function [H, HI, H2] = GABBAEncodedChannelMatrix (h, K) 

% Input: vector of numeric or symbolic entries corresponding to channel estimates 
% Output: GABBA encoded channel matrix H 

% Example: syms hi h2 h3 h4 h5 h6 h7 h8 h9 hlO hll hl2 hl3 hl4 hl5 hl6 
% h = [hi h2 h3 h4 h5 h6 h7 h8 h9 hlO hll hl2 hl3 hl4 hl5 hl6] 

% H = GABBA_EncodedChannelMatrix (h) 

if nargin == 2, 

for k = 1:K % Must be a power of 2 
eval(['syms h' num2str (k) ] ) ; 

end 

h = []; 

for k = 1:K, 

eval(['h= [h; h' num2str(k) '];']); 

end 
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else 

K = length (h) ; % Must be a power of 2 

end 

HI = reshape (h, [1 1 K] ) ; 

H2 = reshape ( [h (K/2 + 1 :end) h(l:K/2)],[l 1 K] ) ; 
while K >= 2; 

for n = l:K/2, 

CI ( : , : , n) = ABBAl (HI ( : , : , 2*n-l) , HI ( : , : , 2*n) ) ; 

C2 ( : , : , n) = ABBA2 (H2 ( : , : , 2*n-l) , H2 ( : , : , 2*n) ) ; 

end 

HI = CI; 
H2 = C2; 
clear CI C2; 
K = K/2; 

end 

K = length (h) ; 

H= [Hl(l:K/2,:) zeros (K/2 , K) ; zeros(K/2,K) H2 (1 :K/2, : ) ] ; 
% Auxiliary functions 

function C = ABBAl (x,y); C = [[ x y] ; [ y -x] ] ; 
function C = ABBA2(x,y); C = [[ x -y] ; [ y x] ] ; 

Appendix III 

Relabeling of GABBA Reduced Encoded Channel Matrices 
A. Matlab Function: Relabels GABBA Reduced Encoded Channel Matrices 
function [Hhatl, Hhat2 ] = RelabelH (HI , H2 , f lag) 

% Input: HI, H2 - minors of GABBA (reduced) encoded channel matrix 

% flag - 1 if minors are from GABBA encoded channel matrix 

% 2 if minors are from GABBA reduced encoded channel matrix 

% Output: relabeled outcome of (sum) product 

g, 
o 

K = length (HI) ; 
switch flag, 
case 1, 

H = conj ( (H1*H1' + H2*H2')/2); 

HI = H (1 :K/2, 1 :K/2) ; 

H2 = H (l+K/2 :K, l+K/2 :K) ; 
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case 2, 

H = HI . ' *H2; 

[pO, pi] = Permutationlndexes ( 1 : K) ; 
HI = H (pO,pO) ; 
H2 = H(pl,pl) ; 

end 

Hhatl = zeros {K/2 ) ; 
Hhat2 = zeros {K/2) ; 
for k = l:K/2, 

eval{['sYms h' num2str (k) ] ) ; 

Hhatl = Hhatl + . . . 

eval ([' double (expand (Hl-Hl (1, ' num2str{k) '))==0).*h' num2str(k) ';'])... 
-eval {[' double (expand (Hl+Hl (1, ' num2str(k) '))==0).*h' num2str (k) ';']); 
Hhat2 = Hhat2 + . . . 

eval ([' double (expand (H2-H1 (1, ' num2str(k) '))==0).*h' num2str(k) ';'])... 
-eval ([' double (expand (H2+H1 (1, ' num2str(k) '))==0).*h' num2str (k) ';']); 

end 

B. Matlab Code: Generator of Relabeled GABBA Reduced Encoded Channel Matrices 

function H = RelabeledGABBAReducedEncodedChannelMatrix (K) 
N = K/2; 
h = []; 

for n = 1:N % Must be a power of 2 
eval(['syms h' num2str (n) ] ) ; 
eval(['h= [h; h' num2str(n) '];']); 

end 

H = reshape (h, [1 1 N] ) ; 
while N >= 2; 

for n = l:N/2, 

C(:, :,n) = ABBA2 (H(:, :,2*n-l) , H ( : , : , 2*n) ) ; 

end 

H = C; clear C; 
N = N/2; 

end 

% Auxiliary functions 

function C = ABBA2(x,y); C = [[ x -y] ; [ y x] ] ; 
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Appendix IV 
Permutation Indexes 

A. Matlab Function: Generator of Permutation Indexes Pq.j^ and Pi-jv 
function [pO, pi] = Permutationlndexes (k) 
% Input: vector of indexes k = 1,2,...K 

% Output: permutation indexes describing the real quasi-orthogonality 
% of GABBA equivalent encoded channel matrices 

% Example: [pO, pi] = Permutationlndexes ( 1 : 32 ) 
N = length (k) ; 
P = k; 

for n = l:ceil (log2 (N) ) -1, 

p = p + floor ( (k-1) / (2^n) ) ; 

end 

p = mod (p, 2 ) ; 

pO = find{p); pi = find(p-l); 

Appendix V 

Example: Orthogonal Decoding of the 4-by-4 GABBA Code 

A. Matlab Script: Explicit Calculations for Orthogonal Decoding of the A-by-A GABBA Code 

% Script: orthogonal soft-decoding of the 4-by-4 GABBA code 
clear all 
K = 4; 

syms si s2 s3 s4 % Symbolic variables declaration 

s = [si s2 s3 s4] .'; % Algebraic symbol vector 

C = GABBAEncoder (s) ; % Algebraic GABBA code (variation of original ABBA) 

H = GABBAEncodedChannelMatrix ( [ ] , K) ; % Algebraic GABBA encoded channel matrix 

r = H* [s; con j (s) ] ; % Algebraic receive vector 

% Linear combinations corresponding to equations 38 and 39 

r_hat(:,l) = collect (collect (expand (H (:, 1 : 2) ' *r + (r' *H ( : , 5 : 6) ) . ' ) , si) , s2) ; 
r_hat(:,2) = collect (collect (expand (H (:, 3 : 4 )' *r + (r' *H ( : , 7 : 8 ) ) . ' ) , s3 ) , s4 ) ; 

% Constructing the GABBA encoded channel matrix algebraically 
HI = H(l:K/2, 1:K); % First minor 

H2 = H((1+K/2):K, (1+K):2*K); % Second minor 
% Algebraic GABBA reduced encoded channel matrix 
HH = simple (conj ( (H1*H1' + H2*H2')/2)); 
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% At this point verify that: 

HH*[sl;s2] == collect (collect (r_hat (:, 1) , si) , s2) % Results ones if true 

HH*[s3;s4] == collect (collect (r_hat (:, 2) , s3) , s4) % Results ones if true 

% Real quasi-orthogonality of GABBA reduced encoded channel matrix 

HHH = HH.'*HH; % Results in a diagonal matrix with identical entries 

% The soft symbol estimates are the result of the linear combinations 

% corresponding to equqations 53 through and 55, normalized by the entries 

% of the diagonal matrix HHH, which are identical 

s_hat = [ ] ; 

s_hat = [s_hat; simple ( (HH .' *r_hat (:, 1 )) /HHH ( 1 ))] ; 
s_hat = [s_hat; simple ( (HH .' *r_hat (:, 2 )) /HHH ( 1 )) ] 

Appendix VI 

Example: Orthogonal Decoding of the 4-by-4 GABBA Code 
A. Matlab Script: Explicit Calculations for the Inspection of Orthogonal Decodability of GABBA Codes 

clear 

% Input the size of GABBA mother encoding matrix 

K = input ('Enter value of K (must be a power of 2): ' ) ; f printf ( ' \n' ) ; 
for n = 1 : log2 (K) -1, 
if n == 1, 

% Generate relabeled algebraic first-order reduced encoded matrix 
HI = RelabeledGABBAReducedEncodedChannelMatrix (K) ; H2 = HI; 

else 

% Generate relabeled algebraic higher-order reduced encoded matrix 
[H1,H2] = RelabelH(Hl,H2,2) ; 

end 

H = HI . ' *H2; 

[pl,p2] = Permutationlndexes ( 1 : K/2 "n) ; 

% Displaying the real-orthogonality of n-th order reduced matrix 

f printf ([' Of f-diagonal Minors of the Real-Product of Permuted GABBA' ... 

' Reduced Encoded Channel Matrix of Order ' num2str(n) ' \n' ] ) 
H(pl,p2) H(p2,pl) 

f printf ([' Type "return" and press ENTER to continue\n' ] ) 

% Waiting for input to continue 

keyboard 

end 
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Appendix VII 
The GABBA Orthogonal Decoder 

A. Matlab Function: GABBA Orthogonal Decoder (Nested Implementation) 

function s_hat = GABBADecoder (rr, hh) 

% Inputs: rr = received signal vector (s) (one comlum per receive antenna) 
% hh = channel estimate vector (s) (one comlum per receive antenna) 

% Output: s_hat = ordered vector of solf orthogonal symbol estimates 
[K,nr] = size(rr); r_hat = zeros (K/2 , 2 ) ; H = zeros (K/2 , K/2 , 2 ) ; 
for a = 1 : nr , 

r = rr ( : , a) ; h = hh ( : , a) ; 

% First Stage 

Hk = GABBAEncodedChannelMatrix (h) ; 

for n = 1:2, 

r_hat(:,n) = r_hat(:,n) + Hk ( : , 1 + (n-l)*K/2 : n*K/2)'*r ... 

+ (r'*Hk(:, 1 + K + (n-1) *K/2 : K + n*K/2)).'; 

end 

HI = Hk (1 :K/2, 1 :K) ; H2 = Hk ( (l+K/2) :K, (1+K) : 2*K) ; 
H(:,:,l) = H(:,:,l) + conj ( (H1*H1' + H2*H2')/2); 

end 

H ( : , : , 2 ) = H ( : , : , 1 ) ; 

clear r rr h hh HI H2 Hk 

% Normalizing Intermediate Quantities 

NormH = norm (H ( : , : , 1 ) ) ; H = H/NormH; r_hat = r_hat/NormH; 

% Higher Stages 

Ps = [[l:K/2].' [(l:K/2) + K/2].']; 
if K > 2, 

[pl,p2] = Permutationlndexes (1 :K/2) ; 
for i = 1 : (log2 (K) -1) , 

r = r_hat ; r_hat = [ ] ; 

Ps_hat = Ps; Ps = [ ] ; 

for n = 1 : 2"i, 

Ps = [Ps Ps_hat (pi (1:K/ (2' (i + 1) ) ) ,n) Ps_hat (p2 ( 1 : K/ (2 " (i+1) ) ) , n) ] ; 
V = H ( : , : ,mod (n, 2) +1) . ' *r ( : , n) ; 

r_hat = [r_hat v (pl(l :K/ (2" (i+1) ) ) ) v (p2 (1 :K/ (2" (i+1) ) ) ) ] ; 

end 



February 1, 2008 



DRAFT 



42 

H_hat = H(:, :,1) .'*H(:, :,2) ; H = []; 

H{:,:,1) = H_hat (pi (1:K/ (2" (i + 1) ) ) ,pl (1:K/ (2" (i + 1) ) ) ) ; 
H{:,:,2) = H_hat (p2 ( 1 : K/ (2 ^ ( i + 1 ) ) ) , p2 ( 1 : K/ (2 ^ ( i + 1 ) ) ) ) ; 

% Normalizing Intermediate Quantities 

NormH = norm (H ( : , : , 1 ) ) ; H = H/NormH; r_hat = r_hat/NormH; 



end 

end 



"6 




for n = 1:K, 




s_hat (Ps { 


n),l) = r_hat (n) /H ( : , : ,mod (n+1, 2) +1) ; 


end 






Appendix VIII 


Exact BER Performance of STBCs in Hoyt, Rice, Rayleigh and Nakagami Fading Channels 


function BER 


= BER_STBC_GeneralFading (esno, Omod, ModType, m, Omega, nr , rho, eta, Naka) 


% Inputs: 




% esno 


= vector of constellation-energy-to-noise-power ratio figures (in dB) 


% Omod 


= number of bits per symbol 


% ModType 


= for QAM, 1 for PSK 


% m 


= vector of "m" parameters (one entry per transmit antenna) 


% Omega 


= vector of "Omega" parameters (one entry per transmit antenna) 


% nr 


= number of receive antennas 


% rho 


= rate of the STBC as defined in Definition 5 


% eta 


= diversity order, as defined in footnote 24 


% Naka 


= leave empty for physically binding, exact, Rice-Hoyt channel model; 


Q, 
O 


or input a dummy value for results with the Nakagami approximation. 


global g phi 


s gamma flag 


flag = 0; 




if nargin == 


9, 


flag = 1; 




end 




nt = length (Omega) ; 


if length (m) 


~= length (Omega) , 


m = m ( 1 ) * 


ones (length (Omega) , 1) ; 


end 
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K = 5000; M = 2"0mod; 
sigma2Vec = (10 . " (-esno/10) ) ; 
if ModType 

% PSK case 
for k = 1:M-1, 

d(k) = 2*abs(k/M - round (k/M)); 
for i = 2 : Omod, 

d(k) = d(k) + 2*abs (k/ (2"i) - round (k/ (2^i) )) ; 

end 

end 

for n = 1 : length (esno) , 

P_minus = []; P_plus = [ ] ; 
for k = 1:M-1, 

% - First integral 

delta = (2*k-l)/M; phi = linspace ( le-50 , pi* ( 1-delta) , K) 
g = ( sin (delta*pi ) " 2 ) /rho; s = -g . / ( sin (phi ) . ' 2 ) ; 
for a = l:nt, 

gamma = Omega (a) / (sigma2Vec (n) ) ; 

MGF_minus (a, : ) = MGF_Explicit (m(a) ) . " (nr*eta) ; 

end 

P_minus (k) = phi (2) *trapz (prod (MGF_minus, 1) ) ; 
% - Second integral 

delta = (2*k+l)/M; phi = linspace (le-50, pi* (1-delta) , K) ; 
g = (sin (delta*pi) "2) /rho; s = -g. / (sin (phi) . "2) ; 
for a = l:nt, 

gamma = Omega (a) / (sigma2Vec (n) ) ; 

MGF_plus (a, : ) = MGF_Explicit (m (a) ) . ' (nr*eta) ; 

end 

P_plus (k) = phi (2) *trapz (prod (MGF_plus, 1) ) ; 

end 

BER(n) = sum (d. * (P_minus - P_plus) ) / (2*pi*0mod) ; 

end 

else 

% QAM case 

phi = linspace (le-50, pi/2, K) ; 
for n = 1 : length (esno) , 
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for k = l:0mod/2, 
I = []; 

for i = : ( (l-l/2"k) *sqrt (M) -1) , 

g= (3* (2*i + l) '2) / (2* (M-1) *rho) ; 
s = -g. / (sin (phi) . "2) ; 
for a = l:nt, 

gamma = Omega (a) / (sigma2Vec (n) ) ; 

MGF (a, : ) = MGF_Explicit (m (a) ) . " (nr*eta) ; 

end 

P(i + l,k) = ( (-1) "floor ( (i*2~ (k - 1) ) / sqrt (M) ) ) * ... 

(2" (k-1) -floor ( (i*2"(k - l))/sqrt(M) + 1/2))* ... 
trapz (prod (MGF, 1) ) ; 

end 

end 

BER(n) = 4*phi (2) *sum (sum (P, 1) , 2) / (sqrt (M) *pi*Omod) ; 

end 

end 

% Moment Generating Functions 

function MGF = MGF_Explicit (m) 
global g phi s gamma flag 
if flag, 

% Nakagami channel 

MGF = (1 - s*gamma/m) . ' (-m) ; 

else 

switch m > 1, 

case % Hoyt channel (reduces to Rayleigh for m = 1) 
q = sqrt((l - 2 * sqrt (m-m" 2 ) ) / ( 2 *m - 1)); 

MGF = l./sqrt(l - 2*s*gamma + ( (2*s*gamma*q) . "2) / ( (l+q"2) "2) ) ; 
case 1 % Rice channel (reduces to Rayleigh for m = 1) 
k = sqrt (m"2 - m) / (m - sqrt (m"2 - m) ) ; 

MGF = ( (1 + k) . / (1 + k - s*gamma) ) . *exp (k*s*gamma. / (1 + k - s*gamma) ) ; 

end 

end 
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Let 



Appendix IX 
Derivation of Equation ( l92t 

1 

Ik = j x)'^'^dx, with fc, X e N. 



By parts we obtain 

h = • (1 - x)^-'^ 
At fc = we have, 

By induction we then find 

h 



k 



K-k+1 



K-k + 1 



fc-(fc-l)---l 



1 {K-k)lkl 



(K - k + 1) ■ {K - k + 2)---K K + 1 (K + iy. 



4-1. 



(95) 



(96) 
(97) 

(98) 

□ 
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